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Abstract 


The  purpose  of  this  research  is  to  conduct  fundamental  investigations  of  tur¬ 
bulent  mixing,  chemical  reaction  and  combustion  processes  in  turbulent,  subsonic 
and  supersonic  flows.  Our  program  is  comprised  of  several  parts: 

a.  an  experimental  effort, 

b.  an  analytical  effort, 

c.  a  computational  effort, 

d.  a  modeling  effort, 

and 

e.  a  diagnostics  development  and  data- acquisition  effort, 

the  latter  as  dictated  by  specific  needs  of  the  experimental  part  of  the  overall  pro¬ 
gram. 

Our  approach  has  been  to  carry  out  a  series  of  detailed  theoretical  and  experi¬ 
mental  studies  of  turbulent  mixing  in  primarily  in  two,  well-defined,  fundamentally 
important  flow  fields:  free  shear  layers  and  axisymmetric  jets. 

To  elucidate  molecular  transport  effects,  experiments  and  theory  concern  them¬ 
selves  with  both  reacting  and  non-reacting  flows  of  liquids  and  gases,  in  fully- 
developed  turbulent  flows,  z.e.,  in  moderate  to  high  Reynolds  number  flows.  The 
computational  studies  are,  at  present,  focused  at  fundamental  issues  pertaining  to 
the  computational  simulation  of  both  compressible  and  incompressible  flows.  Mod¬ 
eling  has  been  focused  on  both  shear  layers  and  turbulent  jets,  with  an  effort  to 
include  the  physics  of  the  molecular  transport  processes,  as  well  as  formulations 
of  models  that  permit  the  full  chemical  kinetics  of  the  combustion  process  to  be 
incorporated.  Our  primary  diagnostic  development  efforts  are  currently  focused 
on  data-acquisition  electronics  to  meet  very  high-speed,  high-volume  data  require¬ 
ments,  the  acquisition  of  single,  or  a  sequence,  of  two-dimensional  images,  and  the 
acquisition  of  data  from  arrays  of  supersonic  flow  sensors.  Progress  has  also  been 
made  in  the  development  of  a  dual-beam  laser  interferometer/correlator  to  measure 
convection  velocities  of  large  scale  structures  in  supersonic  shear  layers  and  in  a  new 
method  to  acquire  velocity  field  data  using  pairs  of  scalar  images  closely  spaced  in 
time. 
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1.  Introduction 

We  have  made  progress  under  the  sponsorship  of  this  Grant,  for  the  period 
ending  14  May  1992,  in  several  areas.  A  review  paper  documenting  much  of  the 
work  at  Caltech  in  shear  layer  mixing  and  combustion  in  the  last  ten  years,  or  so, 
has  appeared  as  Chapter  5  in  the  High-Speed  Flight  Propulsion  Systems ,  Volume 
137  of  the  AIAA  Progress  Series  in  Astronautics  and  Aeronautics  (Dimotakis  1991). 
Copies  are  available  on  request. 

In  our  experimental  effort  in  supersonic  shear  layers,  we  have  increased  the 
high-speed  stream  Mach  number  to  2.5  and  performed  a  first  set  of  exploratory 
investigations  at  this  higher  compressibility  condition.  We  have  also  performed  a 
first  set  of  experiments  in  the  transonic  shear  layer  regime. 

In  our  investigations  of  turbulent  jets,^  two  theses  are  now  available.  An  Exper¬ 
imental  Investigation  of  Chemically-Reacting,  Gas-Phase  Turbulent  Jets  (Gilbrech 
1991),  that  explored  Reynolds  number  effects  on  turbulent  flame  length  and  the 
influence  of  buoyancy  on  turbulent  jet  flames;  and  Mixing  in  High  Schmidt  Num¬ 
ber  Turbulent  Jets  (Miller  1991),  that  examined  Reynolds  number  effects,  Schmidt 
number  effects,  and  the  influence  of  initial  conditions  on  the  turbulent  jet’s  mix¬ 
ing  behavior.  Progress  has  also  been  made  in  understanding  buoyancy  effects  in 
turbulent  jets,  with  an  eye  in  understanding  the  mixing  enhancement  that  we  have 
documented  occurs  under  certain  conditions.  Finally  the  development  of  jet  mixing 
models  that  permit  the  inclusion  of  full  chemical  kinetics  calculations  has  also  pro¬ 
gressed,  with  considerable  improvements  in  the  quantitative  descriptions  of  a  host 
of  phenomena,  including  N0X  production  by  turbulent  jet  flames. 

In  our  analytical/computational  effort,  we  are  developing  new,  reliable  methods 
for  computing  interactions  of  shock  waves  and  other  discontinuous  waves  that  may 
arise  in  multidimensional,  unsteady  compressible  flows. 

In  our  diagnostics  effort,  we  are  continuing  our  progress  with  a  new  generation 
of  high-speed,  high- volume  data  acquisition  systems,  with  the  development  of  a  new 
two-beam  laser  correlator  to  measure  turbulent  structure  convection  velocities,  and 
a  new  method  to  measure  velocity  field  information  from  pairs  of  images. 


+ 


The  investigations  of  turbulent  mixing  and  combustion  in  turbulent  jets  are  co-sponsored  by 
the  Gas  Research  Institute. 
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2.  Mixing  and  combustion  in  supersonic,  turbulent  shear  layers 

Some  preliminary  runs  were  made  during  the  last  year  at  a  high-speed  stream 
Mach  number  of  2.5  and  at  compressible  but  subsonic  flow  conditions.  Both  of 
these  showed  new  and  interesting  phenomena.  However,  on  28  June  1991,  the 
Sierra  Madre  earthquake  caused  a  major  disruption  to  the  Supersonic  Shear  Layer 
(S3L)  combustion  facility.  The  high  pressure  tank  seal  developed  a  leak,  forcing  a 
significant  design  reevaluation  and  repair  effort.  While  that  effort  was  in  progress, 
a  new  instrument  designed  to  measure  large-scale-structure  convection  velocities 
was  being  developed.  These  are  important  for  modeling  chemically  reacting  shear 
layers. 


2.1  Preliminaries  at  Mach  2.5 

The  expansion  of  the  operating  range  of  the  facility  to  a  high-speed  stream 
Mach  number  of  2.5  was  begun  in  May  1991  with  Dr.  Jeffery  Hall.  This  work  was 
left  in  the  preliminary  stage  due  to  several  factors;  the  discovery  of  a  flaw  in  the 
original  design  of  the  S3L  test  section  splitter  plate  and  the  earthquake  damage  to 
the  high  speed  supply  tank.  Despite  these  shortcomings,  significant  new  information 
was  gleaned  from  these  preliminary  runs.  Based  on  previous  experience,  the  growth 
rate  of  these  cases  should  be  very  close  to  runs  at  ideal  conditions;  additionally,  the 
structures  indicated  by  the  travelling  waves  are  a  robust  manifestation  of  supersonic 
convection  velocities  of  turbulent  structures  with  respect  to  the  shear  layer  free 
streams. 

The  new  M\  =  2.5  [N2]  runs  provide  highly  compressible  flows  with  a  different 
density  ratio  from  our  previous  M\  =  1.55  [He]  runs  (Hall  1991).  These  flows  have 
similar  isentropic  convective  Mach  numbers  with  respect  to  each  stream.  Both  runs 
exhibit  travelling  waves  in  the  low-speed  stream,  but  the  wave  speeds  are  differ¬ 
ent,  due,  in  large  part,  to  the  lower  (absolute)  speed  of  the  M\  —  2.5  [N2]  stream 
({7j  =  600m/s),  as  opposed  to  the  M\  =  1.55  [He]  stream  {U\  =  llOOm/s).  The 
normalized  growth  rates  also  are  consistent  (for  the  most  part)  with  the  established 
body  of  data,  at  this  time,  except  for  the  N2/Ar  case.  The  latter,  however,  was 
found  to  have  a  non-zero  pressure  gradient.  Although  the  lowest  compressibility 
case  seems  to  have  a  low  normalized  growth  rate,  this  was  not  unexpected.  By  this 
criterion,  it  is  similar  to  the  low  density  ratio  cases  from  previous  experiments  at 
Mi  ~  1.5. 


FlG.  1  Supersonic  shear  layer  schlieren  image  data:  M\  =  2.5  [NX]  over  L  2  = 
SO  111/s  [Ar]. 

The  .A/]  =  2.5  [He]  runs  showed  some  unexpected  but  not  unprecedented  re¬ 
sults.  These  runs  were  characterized  by  much  higher  compressibility,  i.e...  high 
convective  Mach  numbers,  than  any  previous  runs  in  the  S3L  facility.  The  ex¬ 
pected  travelling  waves  on  the  low-speed  side  were  seen,  and  fast -response  pres¬ 
sure  transducers  were  used  to  measure  their  convection  velocities  using  correlation 
techniques.  Unexpectedly,  the  high  speed  stream  also  showed  signs  of  travelling 
waves.  See  Fig.  2.  The  fast  pressure  transducers  confirmed  this,  giving  a  velocity  of 
t rcj  ~  350 m/s  in  a  stream  with  a  velocity  close  to  U 1  =  1500m/s. 

Two  other  features  of  interest  were  observed  in  the  same  cases.  A  density 
decrease  on  the  upper  edge  of  the  layer  ( c.f.  Figs.  1  and  2)  and  a  normalized  growth 
rate  of  30  —  35f/f  of  the  incompressible  value.  Most  presently  available  data  at  such 
high  compressibility  show  a  growth  rate  only  ~  20(7(  of  the  incompressible  value. 
The  density  decrease  is  of  interest  because  it  implies  that  the  helium  is  heating  up 
from  the  freest  ream  (static)  temperature  value  of  T\  ~  97  Iv  before  mixing  to  any 
significant  extent  with  the  NX  low-speed  stream.  There  is  reason  to  believe  that  the 
three  unusual  factors  are  related. 

The  schlieren  data  picture's  of  these  flows  show  evidence  of  the  original  design 
shortcomings.  Waves  can  be  seen  to  emanate  from  the  supersonic  nozzle  in  both 
cases.  In  the  X>  run.  this  was  caused  by  one  factor.  The  original  design  underesti¬ 
mated  the  bending  loads  on  the  test  section  splitter  plate  tip:  the  support  is  too  far 
back  to  prevent  the  plate  from  elastically  bending  by  about  O.OGOin..  when  exposed 
to  the  y >0 1  ~  17  atm  required  in  the  plenum  for  M\  ~  2.5.  The  deflection,  in  effect, 
results  in  a  two-dimensional  expansion  nozzle  with  a  slightly  incorrect  shape. 


This  problem  has  been  addressed  by  a  redesign  that  should  result  in  a  deflection 


FlG.  2  Supersonic  shear  layer  schlieren  image  data:  M\  ~  2.5  [He]  over  Ui  — 
100  m/s  [N2]. 

that  is  less  than  0.001  in.,  with  plenum  pressures  as  high  as  50  atm.  By  way  of 
comparison,  the  M\  ~  1.5  flows  require  a  plenum  pressure  that  is  less  than  4 atm. 
In  the  helium  high-speed  stream  case,  the  waves  are  stronger  and  the  shear  layer 
can  be  seen  to  bend.  This  is  a  result  of  using  the  same  nozzle,  originally  designed 
for  N2  (7  =  cp/c,.  =  7/5),  with  He  (7  =  5/3),  as  we  did  for  the  A/j  =  1.5  runs. 
At  Mi  =  1.5.  the  two  contours  are  almost  identical.  At  Mi  =  2.5.  however,  the 
differences  are  too  large.  A  new  nozzle  block  has  been  designed  and  fabricated  to 
rectify  this  problem. 

In  previous  experiments  in  the  S  5L  facility,  we  had  found  that  even  strong  wave 
interactions  with  the  shear  layer  produced  relatively  weak  effects  on  the  growth  rate 
of  the  shear  layer.  This  gives  some  confidence  that  the  growth  rate  observed  in  the 
current  flows  is  close  (within  a  few  percent)  to  the  growth  rate  of  a  compressible 
shear  layer  that  should  be  expected  under  ideal  flow  conditions.  The  M\  —  2.5  [N2] 
flows  are  only  slightly  affected;  the  shear  layer  is  still  straight.  The  M 1  =  2.5  [He] 
flows  are  more  affected  by  the  nozzle  shape  mismatch:  the  shear  layer  can  be  seen 
to  bend  on  interaction  with  the  waves  emanating  from  the  nozzle.  Despite  this,  the 
turbulence-generated  travelling  waves  are  a  robust  feature  of  the  flow.  They  initially 
appear  at  the  splitter  plate  and  do  not  change  as  the  stationary  waves  emanating 
from  the  nozzle  interact  with  them.  As  a  consequence,  we  can  expect  that  the 
growth  rate  (Fig.  3)  and  the  travelling  waves  are  not  significantly  influenced  by  the 
complications  resulting  from  nozzle  shape  imperfections. 
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FlG.  3  Compressible  shear  layer  normalized  growth  rate  data:  New  i\/j  =  2.5  and 
subsonic  runs,  as  well  as  previous  data. 

2.2  Compressible  subsonic  shear  layer  flows 

Several  runs  were  made  with  a  transonic,  Ad\  =  0.97  [He]  high-speed  stream, 
over  a  I/2  =  50m/s[N2]  low-speed  stream.  As  can  be  seen  in  Fig.  4,  although  both 
sides  were  subsonic,  weak  travelling  waves  are  discernible  on  the  low-speed  side. 
It  should  be  noted  that,  under  these  conditions,  the  growth  rate  is  so  high  that 
the  shear  layer  attaches  to  the  lower  guidewall  (just  beyond  the  end  of  the  field  of 
view  in  the  picture).  The  contrast  between  the  shear  layer  in  Fig. 4  and  those  in 
Figs.  1  and  2  should  not  escape  unnoticed.  Despite  this  high  growth,  the  normalized 
growth  rate  fits  well  with  other  cases  of  comparable  compressibility.  This  addresses 
a  significant  issue  and  provides  evidence  that  the  whether  the  flow  is  hyperbolic,  or 
not  ( e.g.,  Papamoschou  1989),  is  less  important  than  the  compressibility  of  the  shear 
layer.  The  growth  rates  for  this,  and  a  less  compressible,  case  with  no  travelling 
waves  are  plotted  in  Fig.  3. 
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FlG.  4  Transonic  shear  layer  sclilieren  image  data:  M\  =  0.97  [He],  over  l'<  = 
50  m/s  [No]. 

2.3  Recent  earthquakes  and  required  repairs 

Work  on  both  the  efforts  described  above  was  interrupted  by  a  major  repair 
job  necessitated  by  the  .Tune  1991  Sierra  Madre  earthquake.  After  discovering  the 
leak,  wo  wore  forced  to  stop  experiments.  In  addition  to  considering  simply  re¬ 
placing  the  seal,  other  alternatives  were  investigated.  These  alternatives  included 
remanufneturing  the  tank  to  eliminate  the  large  30 in.  flange.  Unfortunately,  the 
necessity  of  including  the  high  heat  capacity  aluminum  packing  to  prevent  a  sig¬ 
nificant  stagnation  temperature  drop  in  the  high  speed  stream,  in  the  course  of  a 
run.  prevented  the  most  favorable  resolution:  the  heat-treating  temperature  of  the 
welds  (as  required  by  the  high  pressure  vessel  code  for  hydrogen-capable  tanks)  is 
too  close  to  the  melting  temperature  of  aluminum. 

In  the  end.  we  decided  to  replace  the  seal,  with  several  stops  taken  to  minimize 
the  probability  of  future  problems.  Specifically,  the  seal  material  was  changed  to 
one  that  gives  the  closest  coefficient  of  thermal  expansion  to  minimize  thermal 
stresses  when  the  tank  heating  is  activated.  In  addition,  a  ultrasonic  measuring 
device  was  used,  in  retoiquing  the  large  flange  bolts,  to  ensure  that  a  uniform 
and  sufficient  tension  was  placed  on  the  bolts.  This  should  give  the  best  chance  to 
prevent  damage  from  any  future  event  by  giving  us  a  level  of  consistency  which  could 
not  be  obtained  merely  by  tightening  to  a  fixed  torque.  For  better  or  for  worse, 
we  had  the  opportunity  to  test  this.  Following  the  recent  Landers  earthquake  ( 28 
.Tune  1992.  magnitude  7.4.  90  miles  east  of  Pasadena)  and  long  sequence'  of  strong 
aftershocks,  we  loaded  the  large  tank  to  working  presstire  with  helium.  Xo  leaks 
were  found  using  a  sensitive  helium  leak  detector. 
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2.4  Design  of  a  new  instrument 

At  this  time,  travelling  waves  provide  us  with  the  only  means  of  measuring  the 
turbulent  structure(s)  convection  velocity.  Accurate  chemical  kinetic  modelling  of 
the  reacting  shear  layer  depends  significantly  upon  this  velocity  as  it  determines  the 
time  of  flight  and  also  enters  in  the  modeling  of  the  entrainment  ratio.  Specifically, 
a  higher  velocity  gives  a  shorter  time  for  reactions  and  thus  raises  the  required  con¬ 
centrations/temperatures  for  fast  chemistry.  A  relative  velocity  ratio  with  respect 
to  the  two  freestreams  provides  an  asymmetric  entrainment  environment.  There  is 
a  possibility  that  the  wave  velocity  might  not  be  the  only  velocity  in  the  shear  layer; 
if  the  wave  velocity  is  not  the  fluid  velocity,  then  a  model  based  on  the  wave  velocity 
could  be  significantly  wrong.  For  this  reason,  a  new  instrument  was  designed  to 
measure  density  (index  of  refraction)  interface  convection  velocities.  This  should  be 
able  to  measure  both  the  fluid  velocity  and  wave  velocities,  whether  they  are  the 
same  or  different. 
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3.  Mixing  and  combustion  in  turbulent  jets 

The  research  effort  on  turbulent  jet  mixing,  involving  both  the  gas-phase,  chem¬ 
ically  reacting  and  liquid-phase,  non-reacting  jet  investigations,  is  cosponsored  by 
the  Gas  Research  Institute,  GRI  Contract  No.  5087-260-1467. 


3.1  Gas-phase  chemically  reacting  jet  investigations 

The  results  described  below  are  from  the  Ph.D.  thesis  by  Gilbrech  (1991)  that 
was  completed  during  the  last  year  (copies  are  available  on  request).  A  paper 
describing  some  of  the  work  (Gilbrech  &  Dimotakis  1992)  is  included  as  Appendix  B 
of  this  report.  The  experimental  investigations  reported  in  the  thesis  investigate 
a  steady,  effectively  unconfined,  chemically-reacting  turbulent  jet  issuing  into  a 
quiescent  reservoir.  The  Reynolds  number  (based  on  the  nozzle  exit  diameter)  is 
varied  through  pressure,  i.  e.,  density,  with  the  jet  exit  velocity  and  exit  diameter 
held  constant.  The  tiine-averaged  line  integral  of  temperature,  measured  along  the 
transverse  axis  of  the  jet  by  the  wires,  displays  a  logarithmic  dependence  on  x/d* 
within  the  flame  zone,  and  asymptotes  to  a  constant  value  beyond  the  flame  tip.  The 
main  result  of  the  work  is  that  the  flame  length,  as  estimated  from  the  temperature 
measurements,  varies  with  changes  in  Reynolds  number,  suggesting  that  the  mixing 
process  is  not  Reynolds  number  independent  up  to  Re  =  1.5  x  105.  Specifically,  the 
normalized  flame  length  L{/d *  displays  a  linear  dependence  on  the  stoichiometric 
mixture  ratio  <p ,  with  a  slope  that  decreases  from  Re  ~  1.0  x  104  to  Re  ~  2.0  x  104, 
and  then  remains  constant  for  Re  >  2.0  x  104.  Additionally,  the  measurements 
revealed  a  “mixing  virtual  origin,”  defined  as  the  far-field  flame  length  extrapolated 
to  d  =  0.  that  increases  with  increasing  Reynolds  number,  for  Re  <  2.0  x  104  and 
then  decreases  with  increasing  Reynolds  number  for  Re  >  2.0  x  104. 

The  chemically  reacting  jet  experiments  were  performed  at  a  Damkohler  num¬ 
ber  greater  than  5.2,  a  value  foun  I  to  be  the  threshold  to  the  fast  chemistry  regime. 

To  insure  that  the  turbulent  jets  investigated  were  momentum  dominated  to  the 
farthest  measuring  station,  a  separate  set  of  experiments  was  performed.  In  these 
experiments,  the  adiabatic  flame  temperature  AT)  was  raised  to  where  buoyancy 
affected  the  product  thickness  6p  distribution  and  the  logarithmic  dependence  on 
x/d*  of  the  time-averaged  line  integral  of  temperature  no  longer  held.  Under  these 
conditions,  it  was  found  that  molecular  mixing  is  enhanced. 
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3.1.1  Stoichiometric  mixture  ratio  and  Reynolds  number  effects 


All  runs  were  performed  at  pressures  between  1  atm  <  p  <  15  atm.  The  jet  gas 
was  always  F2  in  a  N2  diluent  and  the  reservoir  gas  was  always  NO  in  a  N2  diluent. 


Experiments  at  four  stoichiometric  mixture  ratios,  <j>,  were  performed  at  each 
Reynolds  number  to  produce  four  different  flame  lengths.  In  this  study,  <f>  is  defined 
in  the  usual  manner  as  the  mass  fraction  of  F2  in  N2,  divided  by  the  mass  fraction 
of  NO  in  N2,  normalized  by  the  corresponding  stoichiometric  mass  fraction  ratio, 


i.e., 


_  aw^No) 

“  (rFj/rNO)st 


(i) 


where  the  mass  fraction,  is  defined  for  a  mixture  of  N  species,  as 


£  m 

1=1 


(2) 


and 


m,  =  mass  of  the  ilh  species  . 


The  adiabatic  flame  temperature  rise  AXf  for  each  0  was  computed  using  the 
CHEMIvIN  chemical  kinetics  program  (Kee  et  al.  1980).  The  four  values  of  0 
with  the  corresponding  mass  fractions  and  ATf’s  are  given  in  the  following  table. 


Table  1:  Stoichiometric  mixture  ratios,  mass  fractions,  and 
adiabatic  flame  temperature  rises. 

0 

Yr, 

Ino 

AXf  (K) 

7.0 

0.0127 

0.0029 

7.2 

10.0 

0.0174 

0.0028 

6.9 

14.0 

0.0237 

0.0027 

6.9 

17.9 

0.0300 

0.0027 

6.7 

The  choice  of  AXf  was  determined  from  several  preliminary  experiments.  The 
lowest  value  of  AXf  that  could  be  measured  accurately  and  with  repeatability  was 
sought  to  reduce  both  buoyancy  effects  and  heat  release  effects  on  the  flow  dynamics. 
The  choice  of  AXf  %  7  K.  on  the  basis  of  buoyancy  concerns,  will  be  developed  below. 
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1.4  1.6  1.8  2.0  2.2  2.4  2.6 

log10(x/d*) 

Fig. 5  Normalized  product  thickness  versus  log10(x/d*)  for  the  various 

Figure  5  shows  the  data  for  the  <j>  —  7,  10,  14,  and  18  runs.  The  solid  line  is 
a  compound  curve  fit  with  a  linear  least  squares  fit  in  the  ramp  region  and  a  cubic 
spline  fit  to  the  knee  region.  The  runs  were  fairly  repeatable  as  shown  by  the  good 
agreement  of  the  repeated  runs.  All  of  the  runs  exhibit  a  linear  ramp  region  when 
plotted  in  semi-logarithmic  coordinates  and  a  nearly  constant  asymptotic  level  after 
the  knee  or  flame  tip  region.  The  average  flame  lengths  were  estimated  from  these 
curves  by  determining  the  99%  of  maximum  6p/Lw  point.  These  flame  lengths  are 
plotted  versus  $  for  the  various  Reynolds  numbers  in  Fig.  6.  The  straight  lines  are 
linear  least-squares  fits  to  the  data  of  the  form 

^  =  Ad>  +  B  ,  (3) 

where  .4  and  B  are  functions  of  Re  with  A  being  the  slope  and  B  being  defined 
as  the  “mixing  virtual  origin”.  This  mixing  virtual  origin  is  a  new  observation  and 
corresponds  to  the  abscissa-intercept  in  Fig.  6,  t.  e.,  the  flame  length  extrapolated 
to  4>  =  0.  If  the  flow  could  be  described  by  its  far-field  behavior  at  the  nozzle  exit, 


12 


FlG.  6  Normalized  flame  length  versus  <f>. 

then  B  would  be  the  flame  length  measured  in  the  limit  of  <f>  — ►  0.  The  variation 
of  B  in  Fig.  6  indicates  that  this  mixing  virtual  origin  systematically  decreases 
with  increasing  Reynolds  number  for  Re  >  2.0  x  104.  Some  of  the  variation  in  B 
could  be  attributable  to  the  uncertainty  in  the  curve  fits  used  for  the  flame  length 
determination.  Another  possibility  could  be  a  variation  in  xo,  the  virtual  origin  of 
the  jet,  with  increasing  Re.  However,  even  combined,  these  effects  are  not  large 
enough  to  explain  such  a  large  variation  in  B.  A  nonzero  intercept  in  itself  is  not 
surprising,  if  one  realizes  that  the  potential  core  extends  roughly  six  diameters  from 
the  nozzle  exit  and  that  the  self-similar  region  of  the  jet  typically  begins  around 
x/do  =  20  (e.g.,  Dowling  &  Dimotakis  1990).  Therefore,  fully-developed  turbulent 
mixing  on  a  molecular  scale  requires  some  distance  downstream  of  the  nozzle  exit 
to  develop,  hence  the  streamwise  delay  in  molecular  mixing.  The  fact  that  B  has 
such  a  large  value  and  that  it  varies  with  Re,  however,  is  a  new  observation. 


AT/AT 
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3.1.2  Kinetics 


The  F2/NO  chemistry  in  this  study  can  be  described  by  the  “effective”  reaction: 


2NO  +  F2  -»  2NOF , 


AQ  =  -75 


kcal 

mole  —  F2 


(4) 


-6  -5  -4  -3  -2  -1 


log10(t)  (sec) 

Fig.  7  Chemical  reaction  times  for  some  of  the  runs  described  in  Sec.  3.1.1,  p  = 
1  atm,  T  =  293  K. 


The  reaction  actually  consists  of  two  chain  reactions  of  different  rates  and  heat 
of  reaction.  To  estimate  the  chemical  reaction  time  rc,  these  reactions  were  investi¬ 
gated  using  the  CHEMKIN  chemical  kinetics  software  package  (Kee  et  al.  1980)  with 
a  constant-mass,  constant-pressure  reactor.  Figure  7  shows  a  sample  result  from 
the  program  for  the  mass  fractions  corresponding  to  the  four  4>'s  used  in  Sec.  3.1.1. 
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FlG.  8  Temperature  doubling  to  check  kinetic  rate,  <f>  —  7,  Uq  —  62.4 m/s.  Dam- 
kohler  number  is  given  for  x/d*  =  30. 

The  reaction  time  rc  was  determined  by  drawing  a  tangent  to  the  curve  at 
its  steepest  point  and  using  the  intersection  of  this  tangent  line  with  the  adiabatic 
flame  temperature  as  the  chemical  time.  The  slowest  rate  is  for  (f>  =  7,  p  =  1  atm 
with  rc  =  0.91  msec. 


For  these  experiments,  it  was  required  that  the  time  lag  between  molecular 
mixing  and  chemical  product  formation  be  negligible.  In  this  mixing-limited  con¬ 
dition,  that  we  will  designate  as  the  fast  chemistry  regime,  there  is  a  one-to-one 
correspondence  between  molecular  mixing  and  product  formation.  The  ratio  of  the 
local  large-scale  passage  time  to  the  chemical  reaction  time  is  one  way  of  estimat¬ 
ing  when  the  reacting  flow  is  in  the  fast  chemistry  regime.  This  ratio  of  these  two 
quantities  is  (one  of)  the  Damkohler  number(s),  defined  by 


(5) 
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where  rg  is  the  local  large-scale  passage  time,  estimated  by  dividing  the  local  jet 
diameter  8(x)  by  the  mean  centerline  velocity  Uc  1,  ».e., 

8(x) 

T6  =  HI 

and  rc  is  the  chemical  reaction  time,  computed  as  described  above.  The  local  jet 
diameter  6(x)  for  a  turbulent  jet  has  been  measured  by  many  experimenters  (e.g., 
White  1974,  p.  509;  Dowling  1988)  and  can  be  approximated  by 

S(x)  as  0.44  x  .  (7) 

An  estimate  of  Uci  is  given  for  a  momentum-dominated,  density-matched  jet  by 
Chen  Sz  Rodi  (1980)  as 

^  =  6.2 _ — _  . 

U0  x  -  x0 

The  jet  exit  diameter  d0  in  the  formula  above  should  actually  be  replaced  by  d*, 
the  momentum  diameter  of  the  jet,  but  for  a  properly  designed  nozzle  with  density- 
matched  jet  and  reservoir  gases,  d*  ss  do.  In  the  limit  of  fast  chemistry,  t.e.,  Da  1, 
chemical  product  formation  tracks  molecular  mixing.  The  lowest  Da  of  the  rims  in 
Sec.  3.1.1  occurs  at  x/d*  —  30  for  the  <t>  —  7  run  of  Fig.  7,  with  a  corresponding 
Damkohler  number  of  Da  =  2.8. 

Since  Uo  was  help  constant  for  the  Reynolds  number  experiments  described 
in  Sec.  3.1.1,  and  the  spreading  rate  of  the  jet  is  assumed  to  be  constant  (t.e., 
Reynolds  number  independent),  rg  for  a  given  x/d*  is  the  same  for  all  of  the  current 
experiments  in  Sec.  3.1.1.  The  chemical  time  rc  decreases  with  reactant  number 
density  (t.e.,  pressure).  This  means  that,  for  the  experiments  described  in  Sec.  3.1.1, 
the  numerator  in  Eq.  5  remains  constant  and  the  denominator  is  decreasing  with 
increasing  Reynolds  number.  Therefore,  demonstrating  fast  chemistry  for  the  lowest 
Reynolds  number  runs  guarantees  that  the  higher  Reynolds  number  experiments  of 
Sec.  3.1.1  are  in  the  fast  chemistry  regime. 

One  way  of  determining  if  a  run  is  in  the  fast  chemistry  regime  is  to  repeat  the 
experiment,  doubling  the  reactant  concentrations  of  both  jet  and  reservoir  fluids. 
This  doubles  the  adiabatic  flame  temperature.  The  kinetic  rate  is  proportional  to 
reactant  number  density  and  also  exponentially  dependent  on  the  reaction  tem¬ 
perature.  Therefore,  a  concentration  doubling  increases  the  kinetic  rate  through 
higher  reactant  concentrations  and  by  increasing  the  reaction  temperature.  The 
temperature  rise  for  the  higher  concentration  run  should  double,  if  the  chemical 
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log10(x/d*) 

Fig.  9  Temperature  doubling  to  check  kinetic  rate.  Damkohler  number  is  given  for 
x/d *  =  30. 

product  formation  is  already  mixing-limited.  If  the  chemistry  was  not  fast  before 
the  concentration  doubling,  the  faster  kinetic  rate  after  the  doubling  should  result 
in  more  than  double  the  chemical  product  and,  correspondingly,  more  than  double 
the  temperature  rise. 

A  similar  temperature  doubling  was  performed  for  the  run  in  Sec.  3.1.1  with 
the  lowest  Damkohler  number,  shown  in  Fig.  S.  These  data  do  not  collapse.  Slow 
chemistry  should  only  affect  the  ramp  region  of  the  curve  and  possibly  result  in 
a  longer  flame  length.  The  asymptotic  level,  however,  should  eventually  agree 
at  some  downstream  location  because  all  of  the  F2  should  react  if  one  goes  far 
enough  downstream.  The  asymptotic  level  does  not  agree  for  the  various  ATf’s. 
The  asymptotic  level  nearly  doubles  for  each  temperature  doubling  leading  to  the 
hypothesis  that  some  of  the  F2  was  being  depleted  by  passivation  of  the  nozzle 
plenum  screens  and  honeycomb  during  each  run.  At  low  concentrations,  doubling 
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Fig.  10  Runs  showing  the  lowest  Damkohler  number  that  still  demonstrates  fast 
chemistry.  Da  is  given  for  x/d*  =  30. 

the  absolute  number  of  moles  of  F2,  by  doubling  the  reactant  concentrations,  should 
reduce  the  passivation  effect  in  half.  This  appears  to  be  what  is  happening  in  Fig.  8, 
supporting  the  passivation  hypothesis. 

The  kinetics  and  the  passivation  problems  disappear  for  the  run  with  ATf  = 
14.4  K  in  Fig.  8,  as  demonstrated  by  the  collapse  of  the  data  with  the  higher 
Reynolds  number  run  that  should  be  in  the  fast  chemistry  regime. 

Figure  9  shows  the  result  of  a  similar  experiment  at  p  =  2  atm  with  two  tem¬ 
perature  doublings.  The  collapse  is  good  and  this  condition  appears  to  be  in  the 
fast  chemistry  regime  with  the  ATf  =  7.2  K  run  having  a  Da  —  5.6. 

Figure  10  compares  the  run  at  p  =  1  atm,  <f>  =  14,  to  the  run  at  p  =  15  atm 
and  the  same  4>  and  flame  temperature.  This  is  the  lowest  value  of  the  Damkohler 
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Fig.  11  Runs  showing  effects  of  slow  chemistry  in  the  ramp  region,  <f>  =  10,  ATf  = 
6.9  K.  Damkohler  number  given  for  x/d*  =  30. 

number  in  all  of  the  experiments  which  still  collapses  well  with  an  experiment  at 
a  much  larger  Damkohler  number.  The  lower  limit  of  Damkohler  number  in  these 
experiments,  in  which  the  flow  can  be  considered  to  be  in  the  fast  chemistry  regime, 
is  therefore  Da  ~  5.2. 

In  the  runs  described  thus  far,  the  Damkoler  number  was  changed  by  varying 
the  chemical  time,  but  changing  the  large-scale  passage  time  r#  is  another  way  of 
varying  Da.  Several  runs  at  different  velocities  and  pressures  were  performed  in 
which  the  Damkohler  number  was  low  but  the  temperature  measurement  beyond 
the  flame  tip  achieved  the  correct  level,  showing  them  to  be  free  of  passivation 
problems.  Figure  11  shows  runs  at  various  velocities  and  Damkohler  numbers. 
The  temperatures  for  the  slow  chemistry  conditions  lag  behind  the  fast  chemistry 
runs  but  eventually  achieve  the  same  level  beyond  the  flame  tip.  This  is  how  slow 
chemistry  was  expected  to  manifest  itself.  The  Damkohler  numbers  are  consistent 
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with  the  limiting  value  of  Da  ~  5.2  discussed  above.  The  lower  velocity  runs  show 
the  effects  of  buoyancy  beyond  the  flame  tip.  The  high-velocity  rim  is  repeatable. 


3.1.3  Buoyancy  Effects 


log10(x/d*) 

FlG.  12  Normalized  product  thickness  versus  log10(i/<i*)  for  several  adiabatic  flame 
temperatures,  He  =  4.0  x  104,  <t>  =  10,  U0  =  62.4m/s,  p  =  4atm. 


Buoyancy  is  an  issue  that  was  explicitly  addressed  in  these  experiments.  The 
positively  buoyant  jet  (t.  e.,  a  jet  in  which  the  momentum  flux  increases  with  a*) 
entrains  at  a  higher  rate  than  the  momentum-dominated  jet  ( e.g .,  Becker  iz  Ya- 
mazaki  1978).  Therefore,  the  similarity  law  that  results  in  a  constant  average 
temperature  measurement  beyond  the  flame  tip  will  change. 
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The  temperature  excess  beyond  the  flame  tip  is  a  conserved  scalar  obeying  the 
similarity  law  for  a  momentum-dominated  jet.  Therefore,  the  average  temperature 
measurement  will  decay  beyond  the  flame  tip  at  some  point  when  buoyancy  forces 
become  comparable  to  the  jet  source  momentum  flux. 


For  the  experiments  described  above,  ATf  was  chosen  by  the  requirement  that 
the  jet  remain  momentum-dominated  up  to  the  last  measuring  station  at  x/d0  = 
240  for  the  (p's  investigated.  This  was  done  by  performing  a  set  of  preliminary 
experiments  at  several  <p' s.  The  longest  flame  length  occurs  at  <p  =  17.9.  This  case 
was  examined  by  systematically  increasing  ATf  until  the  temperature  measurement 
beyond  the  flame  tip  was  no  longer  constant,  but  was  observed  to  decrease.  As 
expected,  the  longest  flame  length  is  the  most  sensitive  to  buoyancy  effects  from 
temperature. 


Figure  12  shows  a  set  of  experiments  spanning  adiabatic  flame  temperatures 
ATf  from  3.5  K  to  215  K.  A  Reynolds  number  of  4.0  x  104,  was  chosen  for  this  scan. 
A  <p  of  10  placed  the  flame  tip  in  the  center  of  the  log10(x/d0)  range,  giving  an  equal 
number  of  points  on  each  side  of  the  knee.  The  temperature  measurements  beyond 
the  flame  tip  reach  a  constant  value  for  the  low  heat  release  runs,  indicating  that 
they  are  momentum-dominated.  As  the  flame  temperature  increases,  however,  the 
temperature  measurements  no  longer  attain  a  constant  level  but  begin  to  decrease 
with  increasing  x/d* ,  as  shown  in  Fig.  12,  for  ATf  =  27.6  K. 


These  data  are  presented  here  as  a  first  effort  to  ascertain  that  the  Reynolds 
number  study  was  conducted  in  the  momentum-dominated  flow  regime.  This  series 
of  experiments  showed  an  enhancement  to  molecular  mixing  under  the  influence  of 
buoyancy.  This  enhancement  depends  on  the  relative  effect  of  buoyancy  to  momen¬ 
tum  forces.  Further  study  of  this  behavior  is  currently  in  progress,  as  this  effect 
plays  an  important  role  in  the  context  of  flow  and  mixing  control. 
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3.2  Liquid-phase  turbulent  jet  mixing 

During  the  past  year,  a  thesis  investigating  the  fine  scale  turbulent  structure 
in  a  high  Schmidt  number  (liquid  phase)  jet  was  completed  (Miller  1991).  Close 
to  one  hundred  copies  have  been  distributed  to  people  throughout  the  world.  Lim¬ 
ited  numbers  are  still  available  upon  request.  An  additional  publication  related  to 
this  work  appeared  in  Physics  of  Fluids  A  (Miller  Sz  Dimotakis  1991b)  during  the 
reporting  period  (Appendix  A  of  this  report),  and  a  third  has  been  submitted  for 
publication  (Miller  Sc  Dimotakis  1992).  Seminars,  lectures,  and  talks  highlighting 
specific  aspects  of  this  work  have  been  given  by  Paul  Miller  and  Paul  Dimotakis 
at  ASU,  Caltech,  U.  Chicago,  UCSD,  USC,  at  the  Institute  of  Advanced  Studies 
in  Princeton,  the  DLR  and  Max  Planck  Institute  (Gottingen,  Germany),  LLNL, 
the  Sandia  STAR  meeting,  and  the  November  1991  meeting  of  the  APS  Division  of 
Fluid  Dynamics. 

Results  to  date  include  an  analysis  of  scalar  interfaces  for  fractal  properties 
(which  differed  substantially  from  previous  reports  in  the  literature),  evidence  of 
Taylor  scales  in  the  turbulent  jet,  and  the  absence  of  a  A:-1  region  in  the  small 
scale  region  of  this  high  Schmidt  number  flow  —  all  of  which  have  been  discussed 
in  previous  reports.  Important  new  observations  concerning  the  applicability  of 
fractal  power-law  ideas  to  turbulence  have  grown  out  of  our  work  in  that  area,  and 
a  paper  discussing  the  issue  is  included  as  Appendix  C.  This  paper,  in  conjunction 
with  Miller  &:  Dimotakis  (1991a),  has  changed  how  fractals  are  viewed  in  the  Fluid 
Dynamics  community. 

Additional  findings  during  this  report  period  include  additional  new  evidence 
that  the  jet  pdf’s,  variance,  and  spectra  vary  with  Reynolds  number ,  and  that  these 
changes  with  Re  are  also  Schmidt  number  dependent.  The  long-time  statistics  of 
the  jet  concentration  field  were  found  to  converge  very  slowly,  scalar  spectra  were 
found  not  to  obey  the  classical  k~ 5/3  power-law  prediction,  a  log-normal  behavior 
was  discovered  at  small  scales  in  the  power  spectra  (with  a  starting  point  which 
scales  like  Re3/4),  and  initial  conditions  in  the  jet  plenum  and  near  the  nozzle  were 
seen  to  influence  the  flow’s  mixing  behavior  in  the  far-field. 

In  the  following  sections,  we  will  describe  some  of  these  new  findings  in  greater 
detail.  For  topics  not  covered  here  and  for  a  more  complete  documentation,  the 
reader  is  directed  to  Dimotakis  &:  Miller  (1990),  Miller  (1991),  as  well  as  Miller 
Dimotakis  ( 1991a, b  and  1992). 


3.2.1  Reynolds  and  Schmidt  number  effects 


An  important  result  from  these  experiments  is  that  the  jet  concentration  field 
changes  with  Reynolds  number.  Reynolds  number  for  these  jets  is  defined  as 


(9) 


where  u0  is  the  jet  exit  velocity,  d  is  the  jet  nozzle  exit  diameter,  and  v  is  the 
kinematic  viscosity  of  the  fluid  (water),  and  was  varied  in  these  experiments  by 
altering  uq.  The  jet  concentration  was  measured  on  the  centerline,  and  the  proba¬ 
bility  density  functions,  or  pdf’s,  were  accumulated.  The  current  results  from  the 
high  Schmidt  number  jet,  for  a  wide  range  of  Re,  are  displayed  in  Fig.  13,  along 
with  data  from  gas-phase  jets  (Dowling  198S). 
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FlG.  13  Comparison  of  current  pdf’s  with  gas-phase  results  from  Dowling  (1988). 


As  is  apparent  in  Fig.  13,  the  high  Schmidt  number  pdf’s  become  taller  and 
narrower  with  increasing  Re,  while  the  gas-phase  pdf’s  vary  little.  It  should  be 
noted  that  a  narrower  pdf  implies  a  mixture  with  a  more  uniform  concentration. 
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i.e.,  more  homogeneous  or  well-mixed.  Therefore,  we  conclude  that  the  behavior 
of  the  pdf’s  in  Fig.  13  suggests  that  the  liquid-phase  jet  is  less  well  mixed  than  the 
gas-phase  jet,  but  that  this  difference  decreases  with  increasing  Re.  Alternatively, 
the  gas-phase  composition  is  relatively  insensitive  to  Re,  while  the  liquid-phase  jet 
becomes  better  mixed  with  increasing  Re. 

Another,  related,  measure  of  the  mixedness  of  the  jet  is  the  scalar,  or  concen¬ 
tration,  variance  defined  as 

_  J*cm»x 

var{c}  =  c'2  =  /  (c  — c)2p(c)dc  ,  (10) 

Jo 

( cf .  Miller  Dimotakis  1991b),  where  emax  is  the  highest  concentration  present  and 
p(c)  is  the  probability  density  function  of  the  jet  fluid  concentration  c,  normalized 
such  that 

/  p(c)  dc  =  1  .  (11) 

Jo 

The  variance  is  typically  nondimensionalized  by  the  mean  squared,  c2.  The  resulting 
quantity  is  a  measure  of  the  width  of  the  distribution  p(c/c),  i.e.,  pdf’s  such  as 
in  Fig.  13.  Just  as  a  narrower  pdf  indicates  a  better  mixed  jet,  smaller  variance 
values  imply  more  uniform  composition,  with  a  perfectly  homogeneous,  constant 
concentration  having  a  variance  of  zero.  The  concentration  variance  was  calculated 
from  the  jet  concentration  data  over  a  range  of  Reynolds  numbers.  The  results  are 
compared  with  the  gas-phase  results  in  Fig.  14. 

The  high  Schmidt  number  points  are  fit  with  a  function  chosen  for  its  simplicity, 
ability  to  approximate  the  observed  behavior,  and  capability  to  either  approach  a 
nonzero  asymptotic  value  (or  zero,  if  it  were  indicated).  The  fit  does  suggest  that 
the  points  are  approaching  a  constant  value  with  increasing  Reynolds  number,  and 
the  determined  asymptote,  0.039,  is  in  good  agreement  with  thr  finding  from  the 
analysis  of  previous  data  of  0.04  (Miller  &z  Dimotakis  1991b),  as  discussed  further 
in  Miller  ( 1991 ). 

The  observed  decrease  of  the  concentration  variance  (or,  equivalently,  the  rms) 
in  the  high  Schmidt  number  jet  is  a  substantial  effect.  The  jet  becomes  better  mixed, 
apparently  approaching  an  asymptotic,  high  Reynolds  number  value  of  the  variance 
and  a  similarity  in  the  concentration  pdf’s.  In  a  span  of  1.5  decades  in  Reynolds 
number,  the  variance  is  reduced  by  over  70%.  The  decline  is  not  from  deteriorating 
resolution,  as  discussed  both  in  Miller  &:  Dimotakis  (1991b)  and  Miller  (1991). 
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FlG.  14  Comparison  of  variance  values  with  gas-phase  results  of  Dowling  Sc  Dimo- 
takis  (1990). 

Like  the  pdf’s,  the  gas-phase  variances  measured  by  Dowling  (1988)  change  lit¬ 
tle  over  the  range  of  Reynolds  numbers  investigated,  indicating  a  fundamental  dif¬ 
ference  between  the  order  unity  Schmidt  number  (gas-phase)  and  the  high  Schmidt 
number  (liquid-phase)  jet  mixing.  This  difference  is  also  reflected  in  the  behavior 
of  the  corresponding  power  spectra  for  the  two  cases.  Figure  15  displays  selected 
concentration  power  spectra  from  the  current  work,  with  the  results  of  Dowling 
(1988)  for  comparison. 

Tne  high  Sc  spectra  have  a  considerably  greater  extent  than  the  order  unity 
Sc  spectra  at  the  high  frequencies.  This  behavior  is  expected  in  view  of  a  species 
diffusion  scale  that  is  smaller  than  the  viscous  diffusion  scale,  which  provides  an 
additional  length  scale  in  the  flow  (as  noted  by  Batchelor  1959).  While  no  constant 
—  1  spectral  slope  is  present,  as  discussed  in  previous  reports,  a  log-normal  range  is 
observed  at  frequencies  higher  than  the  viscous  scale,  suggesting  that  the  particular 
assumptions  of  Batchelor’s  theory  which  lead  to  the  —1  slope  require  reexamination 
(c/.  Dimotakis  Sc  Miller  1990).  These  findings  were  presented  at  the  1991  meeting 
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Fig.  15  Comparison  of  power  spectra  with  gas-phase  results  from  Dowling  (19S8). 

of  the  APS  Division  of  Fluid  Dynamics,  and  are  addressed  at  greater  length  in  a 
recently  submitted  paper  (Miller  &:  Dimotakis  1992). 

In  addition  to  the  greater  extent  of  the  high  frequency  portion  of  the  high  Sc 
spectra,  it  is  also  noted  that  the  power  spectral  densities  are  slightly  larger  than 
the  gas-phase  values,  over  the  entire  range  of  frequencies.  A  possible  factor  in 
this  difference  is  that  Dowling’s  jet  nozzle  was  externally  contoured,  and  a  mild 
co- flow  was  provided  to  satisfy  the  entrainment  requirements  of  the  (enclosed)  jet. 
There  was  no  co-flow  in  the  current  work.  If  the  effect  of  the  co-flow  was  to  reduce 
fluctuations  near  the  nozzle,  and  consequently  throughout  the  jet  (in  light  of  findings 
described  below),  it  may  have  resulted  in  the  broadly  lower  spectral  densities. 
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3.2.2  Effects  of  initial  conditions 

It  became  apparent,  during  the  course  of  the  experiments  discussed  above, 
that  something  was  influencing  the  behavior  of  the  scalar  variance.  The  cause  was 
determined  to  be  the  fluid  disturbance  levels  in  both  the  jet  plenum  and  beneath 
the  nozzle  exit.  Additional  runs  were  conducted  specifically  to  address  this  issue. 
One  type  of  run  made  was  termed  a  “quick  run,”  in  which  the  jet  plenum  was  filled 
rapidly  and  the  run  started  within  a  minute.  Eliminating  the  usual  half-hour  period 
wait  for  the  plenum  fluid  to  settle  resulted  in  a  substantially  higher  level  of  fluid 
motion,  or  disturbances,  within  the  plenum  during  the  run.  A  cruciform  originally 
inside  the  plenum  to  damp  such  fluid  disturbances  was  removed,  except  for  one 
comparison  run.  Along  with  the  “quick  runs,”  other  “quiet  runs”  were  made  in 
which  the  plenum  fluid  was  left  to  settle  for  over  half  an  hour,  allowing  the  fluid 
disturbances  to  decay.  Finally,  one  “quiet  run”  was  made  in  which  the  tank  fluid 
immediately  beneath  the  nozzle  was  given  a  vigorous  stir  with  a  hand,  just  before 
initiating  the  run.  This  is  the  “extra  disturbance”  run  shown  in  Fig.  16,  along  with 
the  other  results. 

It  can  be  seen  in  Fig.  16  that,  within  some  scatter  of  the  results,  the  quick 
runs  have  higher  variance  values  than  their  corresponding  quiet  runs.  The  extra 
disturbance  case  also  has  a  higher  variance,  while  the  rim  with  the  cruciform  in 
place  has  a  low  variance  which  agrees  well  with  a  fit  (dashed  line)  to  the  data  previ¬ 
ously  displayed  in  Fig.  14.  These  quantitative  measures  of  the  effect  of  qualitatively 
different  disturbance  levels  lead  us  to  conclude  that  the  jet  variance  increases  with 
increasing  disturbance  level,  both  in  the  plenum  and  in  the  vicinity  of  the  nozzle 
exit. 


►  This  is  a  somewhat  surprising  finding.  It  says  that  adding  perturbations  (in 

the  form  of  disturbances)  to  the  jet  conditions  causes  the  variance  to  increase ,  i.e., 
making  the  jet  less  well  mixed.  It  raises  the  possibility  that  the  natural,  undisturbed 
turbulent  jet  may  provide  the  most  “efficient”  mixing.  These  results  are  for  the 
(  liquid-phase  (high  Schmidt  number)  jet,  and  should  be  viewed  in  light  of  the  earlier 

discussion  of  Sc  effects.  At  this  point,  the  corresponding  gas-phase  behavior  has 
not  been  examined. 
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FIG.  16  Measured  variance  values  for  the  x/d  —  295  runs.  The  shapes  of  the  symbols 
designate  the  run  conditions.  Additionally,  the  filled  symbols  were  selected 
as  representative  of  an  upper  bound  of  the  measured  values,  and  fit  with 
an  expression  of  the  form  A  +  B  Rep.  Parameter  values  determined  from 
the  fit  are:  A  =  0.045,  B  =  19.6,  and  p  =  —0.54. 

3.2.3  Effects  of  plenum/nozzle  geometry 

A  second  factor,  the  plenum/nozzle  geometry,  was  also  found  to  affect  the 
resulting  scalar  variance  values  in  the  high  Schmidt  number  turbulent  jet.  Earlier 
experiments,  described  in  Miller  Sz  Dimotakis  (1991b)  and  referred  to  in  Miller 
(1991)  as  “Phase  I  data,”  also  displayed  a  decrease  in  the  jet  scalar  variance  as  seen 
in  Fig.  14,  but  had  smaller  values  and  approached  their  asymptotic  behavior  sooner 
(at  a  lower  Reynolds  number).  The  Phase  I  data  are  displayed  in  Fig.  17  as  inverted 
triangles.  They  appear  to  lie  on  a  distinct  “branch”,  apart  from  the  other  data. 

We  attribute  the  difference  in  behavior  seen  in  Fig.  17  to  the  only  parame¬ 
ters  that  were  changed  between  the  two  cases:  the  plenum/nozzle  geometry.  The 
Phase  I  experiments  utilized  a  smaller  diameter  plenum  assembly  than  the  following 
(Phase  II)  experiments,  although  the  nozzle  exit  diameter  was  the  same.  Therefore, 
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Fig.  17  Measured  variar^e  values  for  three  major  types  of  runs.  The  inverted  tri¬ 
angles  represent  a  different  plenum/nozzle  geometry.  The  squares  are  the 
data  displayed  in  Fig.  16.  The  other  symbols  axe  the  data  shown  in  Fig.  14. 

the  contraction  ratios  were  quite  different,  with  Phase  I  being  smaller.  In  addition, 
the  details  of  the  nozzle  geometries  differed  between  the  two  cases.  The  Phase  I 
nozzle  consisted  of  a  smooth  radius  at  the  inlet  and  a  short  length  of  0.254  cm 
(0.1  in.)  I.D.  tubing  extending  l.Scm  (0.7in.)  beyond  the  plenum  endplate.  The 
tubing  extension  was  tapered  on  the  outside  to  provide  a  thin  edge  at  the  nozzle 
exit.  The  Phase  II  nozzle  consisted  of  a  circular  contour  leading  to  the  exit  plane  of 
the  nozzle  plate,  with  a  ~  3°  converging  angle  at  the  exit.  With  no  extension  tube, 
the  nozzle  exit  was  flush  with  the  surface  of  the  lucite  nozzle  plate.  Photographs  of 
the  plenum/nozzle  assemblies  may  be  found  in  Miller  (1991). 

The  conclusion  from  these  results  is  that  the  jet  plenum/nozzle  geometry  affects 
the  mixing  occurring  in  these  turbulent  jets,  as  demonstrated  by  the  changes  in  the 
jet  scalar  variance.  At  this  time,  we  have  not  isolated  which  particular  aspects  of 
the  differences  between  the  geometries  were  responsible  for  this  behavior.  However, 
considering  the  above  discussion  of  the  influence  of  disturbances,  both  the  change  in 


contraction  ratios  and  the  variations  in  the  nozzles  are  consistent  with  that  effect. 
The  contraction  ratio  serves  to  amplify  disturbances  in  the  plenum  fluid  as  it  passes 
into  the  nozzle.  Therefore,  it  might  be  expected  that  a  larger  contraction  ratio 
would  exhibit  a  larger  variance,  which  was  indeed  the  case. 

The  small  extension  on  the  Phase  I  nozzle,  in  contrast  to  the  flush-mounted 
hole  of  the  Phase  II  nozzle,  permitted  a  different  (entrained)  flow  in  the  vicinity 
of  the  nozzle  exit.  Such  a  flow  could  be  more  stable  to  small  disturbances  in  the 
vicinity  of  the  nozzle  exit.  Recalling  that  introducing  a  disturbance  beneath  the 
nozzle  exit  increased  the  variance  (Fig.  16),  it  is  likely  that  the  Phase  I  nozzle  is 
less  sensitive  to  disturbances,  and,  therefore,  yields  a  lower  variance.  This  is  also 
consistent  with  the  observed  behavior. 


3.3  Joint  Caltech-Sandia  National  Laboratories  modeling  effort 

As  has  been  described  in  previous  reports  and  papers,  the  Two-Stage  La- 
grangian  model,  under  further  development  in  collaboration  with  Sandia,  has  cor¬ 
rectly  predicted  the  influence  of  the  several  parameters  that  control  combustion 
generated  NOx.  The  absolute  level  of  the  predicted  emission  index,  however,  has 
been  too  high  by  as  much  as  an  order  of  magnitude.  The  source  of  the  error  has 
now  been  identified. 

The  model,  is  in  essence,  a  description  of  the  nature  and  rate  of  molecular 
mixing,  with  the  rate  determined  from  scaling  arguments  and  constants  taken  from 
experiment.  In  the  past,  the  mixing  rate,  because  of  the  lack  of  other  data,  was 
inferred  from  the  dependence  of  the  flame  length  on  the  stoichiometric  ratio.  Fixing 
the  scaling  law  constants  by  this  procedure  obviously  leads  to  the  correct  dependence 
of  the  flame  length  on  the  various  parameters,  but  since  this  length  is  determined 
by  the  last  element  of  nozzle  fluid,  to  reach  the  stoichiometric  ratio,  the  absolute 
value  of  the  mixing  rate  so  determined  is  much  too  low.  Since  the  N0X  emission 
index  depends  critically  upon  residence  time  which  in  turn  depends  inversely  on  the 
mixing  rate,  it  is  now  clear  that  a  low  rate  yields  too  liigh  an  index. 

What  is  instead  the  slowest  rate  (which  determines  the  flame  length)  is  a  mean 
rate.  For  a  cold  jet,  the  data  from  Dowling  &  Dimotakis  (1990)  provide  the  neces¬ 
sary  information.  For  reacting  jets,  the  temperature  profiles  measured  by  Becker  & 
Yamazaki  (1978),  combined  with  an  energy  conservation  analysis,  yield  the  required 
rate. 
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FlG.  18  Variation  of  NOx  emissions  indices  and  molar  N02/N0X  ratios  with  N2 
dilution  of  the  fuel.  Solid  curves  and  triangles  denote  CO/H2  flames;  dashed 
curves  and  circles  denote  CH4  flames.  Radiant  fractions  for  the  undiluted 
CO/H2  and  CH4  flames  are  0.04  and  0.115  respectively,  decreasing  to  0.01S 
and  0.084  at  the  maximum  dilution. 

A.  E.  Lutz,  Sandia,  has  used  the  new  entrainment  constant  to  compute  the  NO 
and  NO2  emission  values,  as  well  as  other  flame  quantities,  for  the  nitrogen  diluted 
CO  —  H2  and  the  CH4  flames  of  S.  R.  Turns,  Pennsylvania  State.  The  model  and 
experimental  NOx  results  are  in  good  agreement  and  a  joint  paper  is  in  preparation. 
Figs.  18,  19  and  20  and  the  following  verbatim  discussion  are  taken  from  this  paper. 


Figure  18  shows  measurements  and  model  predictions  of  NOx  emission  indices 
(EINOx).  For  the  CO/H2  flames,  E1N0X  decreases  with  dilution  as  a  result  of 
both  lower  temperatures  and  shorter  residence  times.  The  influence  of  the  radiant 


O-atom  Mole  Fraction  xIO3  NO  Rate  (mole/cm3 s)  xIO7  Temperature  (K) 
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Fig.  19  Solutions  for  temperature,  NO  production  rate,  and  0  atom  concentration, 
vs.  axial  location  illustrating  chemical  effects  on  NO  production  rates  in 
CO/H2  flames.  Dashed  lines:  equilibrium  solutions  for  instantaneous  (lo¬ 
cal)  mixture  fractions.  The  superequilibrium  O  atom  fraction  is  larger  in 
the  diluted  flame,  because  the  residence  time  is  a  factor  of  two  smaller. 


energy  loss  is  small  in  these  flames.  In  contrast,  radiant  losses  are  important  in 
the  CH4  flames,  as  witnessed  by  the  adiabatic  NO*  predictions  being  twice  the  non 
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Fig.  20  Solutions  for  the  temperature,  N0X  emission  index,  and  integrated  N2  pro¬ 
duction  vs.  axial  location  in  the  undiluted  CH4  flame.  Solid  curves:  ho¬ 
mogeneous  reactor.  Dashed  curves:  flame-sheet  reactor.  The  integral  is 
weighted  by  the  relative  mass  in  the  reactors,  by  /(w/pu)(mi/m)dx.  w  is 
the  mass  of  N2O  produced  per  unit  volume  in  the  ith  reactor,  u  is  the 
average  velocity.  The  reactions  are:  Thermal,  N2  +  0  — >  NO  +  N;  prompt, 
N2  +  CH  — »  HCN  +  N;  and  the  N2O  route,  N2  +  O  +  M  — >  N2O  +  M 
minus  N2O  +  H  — *  N2  +  OH. 


I 


ft 


ft 


p 


p 


p 


» 


p 


p 


p 


33 

adiabatic  values.  Remarkably,  the  CH4  emission  indices  decrease  only  slightly  with 
dilution,  because,  as  noted  above,  the  opposing  effects  of  dilution  leave  the  flame 
temperature  nearly  constant. 

Analysis  of  the  chemical  rate  processes  leading  to  NO  production  shows  that 
the  thermal  (Zeldovich)  mechanism  is  the  only  significant  route  for  NO  formation 
in  the  CO/H2  flames.  The  rate-detex-  ring  step  in  this  mechanism  is  the  reaction 

N2  +  0  ->  NO  +  N  . 

Since  the  N2  mole  fraction  remains  nearljr  constant,  NO  production  in  these  flames 
depends  on  two  factors:  the  O  atom  concentration  history  and  the  temperature 
history. 

Figure  19  illustrates  the  degree  to  which  predicted  0  atom  concentrations  are 
in  excess  of  their  equilibrium  values  (dashed  curves).  The  0  atom  superequilibrium 
ratio  is  greater  for  the  N2  diluted  flame  ([0]/[0]eq  =  6-2  vs.  1.5),  since  residence 
time  in  this  flame  is  less  than  half  in  the  undiluted  flame  (1.7ms  vs.  4ms).  The 
spikes  in  the  kinetic  0  atom  concentrations  occur  when  the  cold  mixture  in  the 
homogeneous  reactor  is  ignited  by  the  hot  stoichiometric  mixture  in  the  flame  sheet 
reactor. 

More  important  to  NO  production  is  the  departure  from  equilibrium  of  the 
flame  temperature.  Peak  temperatures  are  suppressed  approximately  50  K  and 
250  K  for  the  undiluted  and  50%-N2  flames,  respectively  (Fig.  19).  The  temper¬ 
ature  during  NO  production  is  lower  in  the  diluted  flame,  so  the  NO  production 
rate  is  slower  despite  the  higher  0  atom  concentrations.  The  temperature  effect 
overwhelms  the  concentration  effect,  because  the  thermal  NO  reaction  has  a  large 
activation  energy  (close  to  75kcal/gmole),  as  noted  by  Correa  et  al. ,  (1984). 

In  the  CH4  flames,  the  NO  formation  chemistry  is  complicated  by  the  contri¬ 
bution  of  the  prompt  mechanism,  as  initiated  by  the  reaction  N2  +  CH  — ♦  HCN  -f 
N.  The  free  N  atom  easily  forms  NO  by  collision  with  either  <  2  or  OH.  The  HCN 
radical  can  lead  to  either  production  or  consumption  of  NO.  The  prompt  initiation 
step  breaks  apart  30%  of  the  N2  that  is  consumed,  compared  to  60%  via  the  thermal 
mechanism  (Table  I  and  Fig.  20).  The  remaining  N2  consumption  (10%)  occurs  via 
the  N20  route, 

N2  -f-  O  -h  M  — *  N02  -f-  M  , 

which  has  a  rate  comparable  to  the  thermal  initiation  step, 


ft 


N2  +  O  -  NO  +  N  , 


Table  I 

Contributions  of  various  kinetic  mechanisms 
to  total  NO  production  in  undiluted  CH4  flames 


N2  consumption  (%  of  total  in  both  reactors): 


Path 

Homogeneous 

Flame-sheet 

Thermal 

15 

45 

N20  route 

4 

6 

Prompt 

0 

30 

NO  production  (%  of  total  in  both  reactors): 


Path 

Homogeneous 

Flame-sheet 

Thermal(a) 

21 

100 

Prompt  (b) 

0 

-18 

N02  conversion 

-12 

4 

HNO  loop 

-3 

8 

(a) N2  +  0  -  NO  +  N; 

N  +  02  — >  NO  +  O; 

N  +  OH  -  NO  +  H 

(b) CN  +  0  -»  C  +  NO; 

HCN  +  0  ->  CH  +  NO; 

HCNO  +  H  -  CH2  +  NO; 

HCNO  +  CO  — ►  HCCO  +  NO; 

C02  +  N  — ►  NO  +  CO; 

NCO  +  O  — ►  NO  +  CO 

except  that  most  of  the  N20  returns  to  N2  by  collision  with  H  atoms.  Of  the  total 
N2  consumed  in  the  simulation,  80%  occurs  in  the  flame-sheet  reactor,  which  is  not 
surprising,  because  the  mixture  here  is  hot  during  the  whole  flame  length.  Similarly, 
94%  of  the  total  NO  production  comes  in  the  flame-sheet  reactor.  The  contribution 
of  the  homogeneous  reactor  occurs  via  the  thermal  mechanism  when  the  reactor 
ignites  at  x/d  =  80  (Fig.  20). 

We  exercised  the  model  without  the  prompt  reactions  included  in  the  mecha¬ 
nism  and  found  that  the  computed  EINOx  was  30%  lower  than  predicted  by  the 
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full  mechanism.  While  this  agrees  with  the  N2  consumption  associated  with  the 
prompt  initiation  step  above,  we  caution  that  superposition  does  not  strictly  apply 
to  the  nonlinear  solution. 

The  evolution  of  the  N0X  emission  index  in  the  simulation  (Fig.  20)  indicates 
that  prior  to  the  ignition  of  the  homogeneous  reactor,  NOx  produced  in  the  flame- 
sheet  is  transported  to  the  homogeneous  reactor.  When  the  buoyant  switch  accel¬ 
erates  the  entrainment  of  cold  surroundings  at  x/d  —  70,  a  period  of  transient  NO 
conversion  to  HNO  occurs  in  the  flame-sheet  (dashed  curve)  via  a  three-body  reac¬ 
tion  with  H  atoms.  Shortly  thereafter,  when  the  homogeneous  reactor  ignites,  NO 
production  proceeds  again,  until  the  fuel  is  exhausted  and  the  temperature  begins 
to  decay.  The  final  emission  level  comes  as  a  result  of  the  fluid  exchange  between 
reactors  as  they  continue  to  entrain  cold  air. 

Next  we  examine  the  NO2  fraction  of  the  total  NOx.  The  NO2  fraction  increases 
as  the  CO/H2  flames  are  increasingly  diluted,  while  the  NO2  fractions  are  unaffected 
by  dilution  in  the  CH4  flame  (Fig.  19).  The  model  predictions  and  experiments  are 
in  excellent  agreement,  especially  considering  the  small  absolute  quantities  of  NO2 
produced.  Although  the  NO2  fractions  increase  in  the  CO/H2  flames,  NO2  emission 
indices  (EIN02)  actually  decrease  with  N2  dilution  (0.137  —  0.061  g/kg  predicted 
for  0-50%  N2).  The  model  predicts  that  the  NO2  is  formed  after  NO  formation  is 
complete,  well  downstream  of  the  maximum  temperature.  It  is  also  interesting  to 
note  that  the  peak  N02  production  occurs  in  all  of  the  flames  in  the  temperature 
range  of  1100  —  HOOK. 

In  conclusion,  predictions  from  a  turbulent  jet  flame  model  that  employs  a 
relatively  simple  model  of  fluid  mixing  but  detailed  chemical  kinetics  were  compared 
with  experimental  measurements  from  ^-diluted  CO/H2  and  CH4  flames.  The 
model  predictions  of  NO  and  NO2  emission  indices  are  in  good  qualitative  and 
reasonable  quantitative  agreement  with  measurements  from  both  flames.  Inclusion 
of  the  radiant  heat  losses  is  necessary  to  predict  NOx  levels  and  trends  with  N2 
dilution  from  the  CH4  flames.  The  thermal  mechanism  is  the  only  significant  rout 
for  NO  formation  in  the  homogeneous  reactor  model  of  CO/H2  flames.  It  is  also 
the  dominant  route  in  the  CH4  flames,  except  for  a  30%  contribution  of  the  prompt 
mechanism,  which  appear  in  the  flame-sheet  reactor  only. 

The  different  channels  of  NOx  formation  reasonably  well  predicted  by  a  detailed 
chemical  kinetics  model  contribute  to  a  better  understanding  and  controlling  NOx 
emissions  from  combustion  sources.  The  model  is  now  in  good  agreement  with 
experimental  results,  after  the  mixing  rate  was  determined  to  be  too  low,  originally. 
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4.  Analytical  and  computational  effort 

A  new  approach  for  computing  multidimensional  flows  of  an  inviscid  gas  has 
been  developed  as  part  of  our  on-going  computational  effort.  It  may  be  considered 
an  extension  of  the  method  of  characteristics  to  multidimensional,  unsteady  gas 
dynamics. 

The  idea  is  to  use  the  knowledge  of  the  one- dimensional  characteristic  problem 
for  gasdynamics  to  compute  genuinely  multidimensional  flows  in  a  mathematically 
consistent  way.  A  family  of  spacetime  manifolds  is  found  on  which  an  “equivalent” 
1-D  problem  holds.  These  manifolds  are  referred  to  as  Riemann  Invariant  Manifolds 
and  play  the  role  that  the  characteristic  manifolds  play  in  the  method  of  character¬ 
istics  for  1-D  unsteady  or  2-D  steady  gas  dynamics.  Their  geometry  depends  on  the 
local  spatial  gradients  of  the  flow  and  they  provide  locally  a  convenient  system  of 
coordinate  surfaces  for  spacetime.  In  the  case  of  zero  entropy  gradients,  functions 
analogous  to  the  Riemann  invariants  of  1-D  gasdynamics  can  be  introduced.  These 
generalized  Riemann  Invariants  are  constant  on  the  Riemann  Invariant  Manifolds. 
The  equations  of  motion  are  integrable  on  these  manifolds  and  the  problem  of  com¬ 
puting  the  solution  becomes  that  of  determining  the  geometry  of  these  manifolds 
locally  in  spacetime.  This  theory  allows  one  to  look  at  such  flows  in  a  very  useful 
way,  reminiscent  of  the  method  of  characteristics.  A  whole  family  of  numerical 
schemes  can  be  devised  using  these  ideas. 

The  progress  made  in  recent  years  in  the  development  of  numerical  schemes 
for  computing  such  flows  is  based  on  the  various  higher  order  extensions  of  the 
original  Godunov  (1959)  scheme.  These  schemes  allow  high  order  accuracy  without 
spurious  oscillations  near  the  flow  discontinuities.  The  main  idea  of  the  Godunov 
scheme  and  its  extensions  is  to  use  the  knowledge  from  the  theory  of  characteristics, 
locally  in  each  computational  cell,  to  compute  the  various  flux  terms.  By  doing 
this,  the  local  characteristic  wave  patterns  are  accounted  for,  accurately.  Since 
discontinuities  are  likely  to  be  present,  the  characteristic  problem  is  generalized  to 
allow  for  their  presence.  This  constitutes  the  well  known  Riemann  problem.  Two 
arbitrary  constant  states,  separated  by  a  discontinuity  in  space,  are  given  as  the 
initial  condition.  The  solution  to  the  Riemann  problem  is  used  as  the  building  block 
for  numerical  methods  of  this  type.  It  provides  the  information  needed  to  compute 
the  fluxes  at  the  interfaces  of  the  computational  cells. 

Although  the  application  of  this  idea  is  straightforward  for  the  one-dimensional 
case,  the  extension  to  multidimensional  flows  is  not.  Practically,  these  schemes  are 


37 


extended  to  more  dimensions  by  treating  the  additional  spatial  dimensions  sepa¬ 
rately.  The  computation  is  not  truly  multidimensional,  but  rather  a  number  of 
separate  one-dimensional  computations  in  each  spatial  direction.  Even  schemes 
that  appear  to  be  multidimensional  (unsplit),  actually  use  the  local  directions  of 
the  grid  to  set  up  and  solve  a  1-D  Riemann  problem,  thus  taking  advantage  of  the 
1-D  results.  The  only  real  effort  to  create  a  multidimensional  scheme  are  those  by 
Roe  (1986),  Hirsch  Sz  Lacor  (1989),  and  Deconick  et  al.  (1986). 

We  have  used  the  Riemann  Invariant  Manifold  theory  to  develop  a  novel 
scheme.  This  scheme  is  a  standard  high-order  shock-capturing  Godunov  scheme 
which  has  been  modified  to  take  into  account  the  multidimensional  character  of  the 
flow  in  a  mathematically  consistent  way.  No  flux-splitting  or  other  multidimensional 
“fixes”,  that  are  usually  employed,  are  needed.  This  multidimensional  correction 
can  be  directly  extended  to  three-dimensional  flows.  We  are  currently  testing  this 
scheme  on  a  variety  of  test  problems. 

A  first  paper,  documenting  our  earlier,  1-D  unsteady  gasdynamic  flow  simula¬ 
tion  method,  has  been  accepted  for  publication  (Lappas  et  al.  1992). 
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5.  Diagnostics  and  instrumentation  development 

Our  diagnostic  and  instrumentation  development  is  proceeding  very  well.  A 
new  instrument  to  measure  the  turbulent  structure(s)  convection  velocity  was  pre¬ 
viously  mentioned  in  Sec.  2.4,  and  will  not  be  discussed  further  here. 

As  part  of  a  major  effort,  we  have  developed  the  capability  to  acquire  data  at 
very  high  rates,  beyond  those  achievable  with  commercial  devices.  This  will  permit 
us  to  acquire  sequences  of  high  resolution  images,  using  CCD  cameras.  Additionally, 
we  are  proceeding  with  the  design  of  electronics  to  support  a  custom  CCD  array, 
designed  in  collaboration  with  JPL,  which  has  been  dubbed  the  “MACH  CCD.”  The 
MACH  CCD  will  permit  us  to  record  two  high-resolution  images  within  a  fraction 
of  a  microsecond  of  each  other.  We  will  then  be  able  to  image  supersonic  flows,  and, 
utilizing  image  correlation  velocimetry  techniques,  determine  the  two-dimensional 
velocity  field  in  the  plane  of  the  images. 


5.1  High-speed  data  acquisition 

We  have  completed  the  first  two  A/D  converter  boards  for  high  speed  and  high 
resolution  image  acquisition.  The  first  board  has  two  12  bit  15  megahertz  A/D 
converters  and  can  digitize  two  analog  inputs  at  15  MHz  or  one  input  at  30  MHz. 
The  second  board  has  two  12  bit  20  megahertz  A/D  converters  and  can  digitize 
two  inputs  at  20  MHz  or  one  input  at  40  MHz.  Both  boards  have  32  megabytes  of 
local  buffer  memory  and  plug  into  a  standard  VXI  backplane.  A  VMEbus  computer 
running  the  OS-9  real  time  operating  system  controls  the  A/D  converters  through 
a  VME  to  VXI  bus  converter.  The  OS-9  system  can  transfer  the  acquired  images 
to  the  VAXcluster  through  the  Ethernet. 


The  A/D  subsystem  has  been  used  with  a  Hitachi  high  resolution  RGB  camera. 
The  camera  has  three  CCD’s  (one  per  color)  with  a  resolution  of  560H  by  480V 
pixels  A  phase  locked  loop  has  been  built  for  the  Hitachi  camera  to  reconstruct 
the  pixel  clock  and  synchronize  the  A/D  converter  boards.  We  can  digitize  up  to 
56  consecutive  frames  when  digitizing  2  of  the  three  colors  or  up  to  28  consecutive 
frames  digitizing  all  3  colors. 
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We  have  also  used  a  high  resolution  black  and  white  camera  manufactured  by 
Texas  Instruments  (resolution  of  1134H  by  480V  pixels).  This  camera  has  a  pixel 
clock  output,  so  a  small  board  was  constructed  to  buffer  the  pixel  clock,  horizontal 
sync,  and  vertical  sync  for  the  A/D  converter  boards.  One  A/D  converter  board  can 
digitize  up  to  28  consecutive  frames  (this  can  be  extended  to  56  frames  by  chaining 
a  second  board  after  the  first). 

Additionally,  we  are  designing  a  controller  board  that  will  be  able  to  control 
the  MACH  dual  image  CCD  camera  being  built  by  JPL.  The  controller  board  will 
have  two  microsequencers.  The  first  (pixel)  microsequencer  will  generate  the  com¬ 
plex  waveforms  needed  by  the  MACH  CCD  camera  and  control  the  A/D  converter 
boards.  The  pixel  microsequencer  can  be  programmed  to  digitize  a  subset  of  the 
image  (region  of  interest)  in  order  to  conserve  memory.  The  second  (real  time) 
microsequencer  is  used  to  generate  the  timing  waveforms  for  the  YAG  lasers,  the 
camera  shutter,  and  other  laboratory  devices.  Both  microsequencers  can  be  trig¬ 
gered  from  an  external  control  signal  or  from  the  other  microsequencer. 

For  instance,  the  real  time  microsequencer  can  wait  for  an  external  trigger  signal 
(e.g.  the  start  of  a  run),  then  fire  the  first  YAG  laser,  tell  the  pixel  microsequencer 
to  store  the  first  image,  then  fire  the  second  YAG  laser,  and  finally  tell  the  pixel 
microsequencer  to  start  digitizing  the  two  images.  The  MACH  CCD  should  be  able 
to  acquire  two  images  as  short  as  200  nanoseconds  apart. 


5.2  Image  correlation  velocimetry 


We  are  developing  the  capability  to  use  the  correlation  of  two  scalar  images 
for  the  purpose  of  motion  detection.  In  particular,  we  can  map  the  deformation 
of  a  flow  field  over  time.  Quantities  like  the  velocity,  vorticity,  and  strain  rate  are 
estimated  directly  from  a  pair  of  scalar  images.  Initial  investigations  designed  to 
tune  this  process  have  examined  the  flow  past  a  circular  cylinder. 


To  measure  the  movement  of  a  scalar  image  over  time,  consider  two  images  or 
scalar  fields  (e.g.,  dye,  temperature,  or  density),  spaced  in  time.  It  is  clear  that 
there  is  a  mapping  taking  the  first  image  to  the  second.  One  such  mapping  is 
simply  the  motion  of  the  fluid  together  with  the  equations  of  motion  of  the  scalar. 
Integrating  the  equation  of  motion  of  a  scalar,  over  time,  we  may  write: 


2(l  (*(x))  -  lto  (x)  -  D 


V2J(x(<))  dt  =  0  , 


(12) 


where  £(x)  is  some  (vector)  coordinate  transformation  mapping  the  image  taken  at 
time  to,  Zto  (x),  to  the  image  at  time  t\,Xix  (x).  The  unknown,  and  highly  nonlinear 
vector  transformation  £(x),  is  then  calculated  by  solving  Eq.  12.  A  more  reliable 
method  is  to  minimize  Eq.  12  over  some  volume  V,  e.g.,  in  a  least  squares  sense, 


min 

€o 


/  Ztl  (* (x))  -  Jl0  (x)  -  V  V2Z  (x(0)  dt 
Jv  L  Jto 


2 

dx 


(13) 


Now  consider  that  £  (x)  is  just  a  displacement  of  some  point  x  to  another  $  (x). 
Since  £  (x)  is  a  complicated  nonlinear  function  we  linearize  (locally  expand)  in  a 
neighborhood  Af  about  a  point  x,  so, 


£(x)  =  £(x)  +  (x-x)'  V£(x)  +  ...  (14) 

Clearly  £  (x)  is  the  displacement  (velocity)  of  the  center  of  „V,  and  (x  —  x)  •  V£ (x)  is 
the  local  deformation  (rotation,  strain,  dilatation,  etc.),  that  Af  experienced  when 
moving  from  x  to  £  (x). 


Note  that  the  state  of  the  ait  in  image  velocimetry  does  no  better  than  the 
first  order  displacement,  £(x),  of  A1".  Using  Eqns.  13  and  14  over  a  neighborhood 
AT  yields, 

min  /  |ztl  (£  (x)  +  (x  —  x)  •  V£  (x)  +  ...)  —  Zt  (x)  — 

£(x),v£(x),...  Jj4  L 

rt\ 

V  /  V2 
Jto 

Note  that  there  is  little  justification  in  going  beyond  the  second  order  term  in  Eq.  15. 
This  is  because  of  the  rapid  increase  in  the  number  of  parameters  that  occurs  when 
going  to  higher  order  in  the  minimization.  Trying  to  fit  too  many  parameters  would 
be  time-consuming,  and,  in  the  presence  of  noise,  meaningless. 

Our  experience  suggests  that  a  second  order  expansion  of  the  displacement 
field  will  yield  better  velocity  estimates  than  results  using  a  traditional  first  or¬ 
der  approximation.  Keeping  the  second  order  terms  in  the  minimization  (Eq.  15) 
increases  the  reliability  of  the  first  order  term,  particularly  in  regions  where  the 
flow  contains  large  velocity  gradients  or  rotational  motions.  In  addition,  the  local 
vorticity,  strain,  and  divergence  can  be  determined  directly  from  the  minimization 
parameters,  obviating  the  large  errors  inherent  in  using  only  the  first  order  estimate 
of  the  displacement  and  then  differentiating  to  obtain  these  higher  order  derivatives 
of  the  flow  field. 


X{x.{t))dt 


dV 


(15) 


The  high-speed  data  acquisition  is  part  of  the  work  by  D.  Lang,  and  the  scalar 
correlation  velocimetry  is  being  carried  out  by  P.  Tokumaru,  both  in  collaboration 
with  P.  Dimotakis. 
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The  scalar  rms  fluctuations  in  a  turbulent  jet  were  investigated  experimentally,  using  high- 
resolution,  laser-induced  fluorescence  techniques.  The  experiments  were  conducted  in  a  high 
Schmidt  number  fluid  (water),  on  the  jet  centerline,  over  a  jet  Reynolds  number  range  of 
3000<Re<24  000.  It  was  found  that  the  normalized  scalar  rms  fluctuations  c'/  c  decrease  with 
increasing  flow  Reynolds  number,  at  least  for  the  range  of  Reynolds  numbers  investigated. 

Since  c' /  c  is  a  measure  of  the  inhomogeneity  of  the  scalar  field,  this  implies  that  high  Schmidt 
number  turbulent  jets  become  more  homogeneous,  or  better  mixed,  with  increasing  Re.  These 
findings  need  to  be  assessed  in  the  context  of  the  documented  Reynolds  number  independence 
of  flame  length  for  Re  >  3000  or  6500. 


I.  INTRODUCTION 

The  correct  description  of  the  probability  density  func¬ 
tion  ( pdf)  of  the  fluctuations  of  a  conserved  passive  scalar  in 
turbulent  flow  has  been  a  matter  of  theoretical  and  engineer¬ 
ing  interest  for  some  time.  Knowing  the  pdf  of  the  conserved 
scalar  allows  the  computation  of  the  amount  of  chemical 
product  that  would  be  formed  in  a  chemical  reaction  in 
which  the  scalar  is  transported  as  a  passive  quantity.' 

Experimental  evidence  and  theoretical  arguments  in  the 
last  few  years  suggest  that  the  turbulent  mixing  process  is 
sensitive  to  the  value  of  the  molecular  diffusion  coefficients, 
even  at  Reynolds  numbers  typically  regarded  as  high.  We 
recognize  that  this  remains  a  controversial  proposal,  at  this 
writing,  especially  in  the  case  of  turbulent  jet  mixing.-1  We 
note,  in  particular,  that  accepting  it  would  imply  that  the 
conserved  scalar  pdf,  as  well  as  other  statistics  of  the  mixing 
process,  should  be  expected  to  be  functions  of  the  flow  Reyn¬ 
olds  number  Re  =  Ub/v,  where  U  is  the  local  flow  velocity, 
<5  is  the  local  extent  of  the  turbulent  flow  and  v  denotes  the 
(molecular)  kinematic  viscosity,  and  the  Schmidt  number 
Sc  =  v/fJ ,  where  fJ  is  the  scalar  species  (molecular)  diffu- 
sivity.  There  is  no  consensus  at  this  time  as  to  how  such 
effects  should  be  described,  let  alone  predicted,  with  the  as¬ 
sociated  nonlinear  dynamics  essentially  involving  terms  in 
the  equations  of  motion  with  factors  that  vanish,  multiplied 
by  factors  that  become  unbounded  in  the  limit  of  Re  —  oc ,  or 
Sc— •  oc . 

Nevertheless,  the  qualitative  dependence  on  Sc  might  be 
argued  on  simple  grounds,  namely,  that  decreasing  the  spe¬ 
cies  diffusivity,  i.e.,  increasing  the  Schmidt  number,  keeping 
all  other  flow  and  fluid  quantities  fixed,  can  only  decrease 
the  rate  of  (molecular)  mixing  and  the  amount  of  mixed 
fluid.  This  is  corroborated  by  a  comparison  of  the  results  of 
experiments  in  gas  phase  shear  layers'  (Sc  —  1 )  versus  those 
in  liquid  phase  shear  layers4  ( Sc  —  10' )  at  comparable  Reyn¬ 
olds  numbers. 

An  a  priori  argument  for  the  dependence  on  Reynolds 
number  is  not  as  straightforward.  Reliable  experimental 


data  at  large  Reynolds  numbers  are  difficult  to  obtain.  Ex¬ 
periments  at  high  Reynolds  numbers  must  rely  on  chemical¬ 
ly  reacting  flows  to  infer  the  degree  of  molecular  mixing.  In 
the  case  of  gas  phase  shear  layers,  the  experimental  evidence 
from  a  single  experiment  is  that  the  amount  of  mixing  de¬ 
creases  as  the  Reynolds  number  increases,  albeit  slowly.5 
The  evidence  shows  a  Reynolds  number  dependence  in  liq¬ 
uid  phase  shear  layers  that  appears  to  be  weaker,  but  this 
conclusion  is  based  on  scanter  data.4 

In  turbulent  jets,  the  behavior  is  even  less  clear.  Experi¬ 
mental  “flame  length”  data  in  liquid  jets,*’  as  well  as  other 
liquid  and  gas  phase  turbulent  jet  flame  length  data  reviewed 
in  Fig.  1  ofDahm  and  Dimotakis,7  suggest  that  this  quantity 
may  well  be  independent  of  both  Schmidt  number  and  Reyn¬ 
olds  number,  at  least  for  Reynolds  numbers  greater  than 
about  3  X  I  O  ',  or  so,  as  had  been  reported  in  the  earlier  Wed¬ 
dell"  (see  also  Hottel,  Ref.  9)  findings.  The  scatter  in  these 
data  is  large  enough,  however,  to  have  masked  a  weak  depen¬ 
dence  on  Schmidt  number  or  Reynolds  number,  not  to  speak 
of  the  likely  contamination  of  the  results  by  buoyancy  effects 
in  the  case  of  high  heat  release,  gas  phase  jet  flames.  As  a  case 
in  point,  the  results  of  Kristmanson  and  Danckwerts"’  sug¬ 
gest  that  Reynolds  number  independence  of  the  flame  length 
in  a  liquid  jet  does  not  occur  until  a  Re  greater  than  6500,  or 
so. 

More  recently,  direct  measurements  of  the  jet  fluid  con¬ 
centration  pdf  in  gas  phase  turbulent  jets  by  Dowling  and 
Dimotakis"  suggest  a  Reynolds  number  independence,  at 
least  in  the  Reynolds  number  range  from  5.0x10'  to 
4.0  X  104.  These  authors  also  compare  their  estimates  of  the 
jet  fluid  concentration  pdf  on  the  jet  axis  with  those  mea¬ 
sured  by  other  investigators,1’1'’  for  Schmidt  (Prandtl) 
numbers  ranging  from  0.7  to  10'  and  Reynolds  numbers 
ranging  from  5x10'  to  5.2 Xl04  (Dowling  and  Dimota¬ 
kis,"  Fig.  28).  While  the  reported  pdf’s  may  not  collapse 
particularly  well,  it  is  difficult  to  assess  whether  the  differ¬ 
ences  are  attributable  to  Schmidt  number  or  Reynolds  num¬ 
ber  effects,  or  to  experimental  difficulties. 

Our  own  investigations  of  the  behavior  of  the  turbulent 
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jet  fluid  concentration  fluctuations  were  undertaken  as  part 
of  a  larger  effort  to  address  the  role  of  Reynolds  number  as  a 
parameter  in  the  dynamics  of  turbulent  jet  mixing.  The  ex¬ 
periments  reported  here  were  performed  in  liquid  phase  tur¬ 
bulent  jets  (Sc~103),  in  a  Reynolds  number  range  from 
3.0  X  103  to  2.4  x  104.  Our  measurements  permitted  direct 
estimates  of  the  pdf,  rms,  as  well  as  the  spectrum  of  the  con¬ 
served  scalar  fluctuations  to  be  computed. 


II.  APPARATUS  AND  PROCEDURE 

While  the  experimental  setup  has  been  documented  pre¬ 
viously,6  17  the  following  short  description  is  included  for 
completeness.  The  experiments  were  carried  out  utilizing 
the  facility  shown  in  Fig.  1.  A  large  rectangular  water  tank, 
about  1  m  square  and  2  m  high,  acted  as  the  (discharge) 
reservoir.  Large  windows  on  all  sides  of  the  tank  provided 
optical  access.  To  establish  the  flow,  the  jet  plenum  was  filled 
with  water  tagged  by  a  fluorescent  laser  dye  (sodium 
fluorescein)  and  air  was  sonically  metered  to  drive  the  jet 
fluid  at  constant  velocity  through  a  2.54  mm  (0.1  in.)  diam 
nozzle  at  the  base  of  the  plenum.  The  beam  from  an  argon- 


ion  laser  was  passed  through  appropriate  optics  and  aligned 
radially  through  the  centerline  of  the  jet.  A  beam  stop  pre¬ 
vented  reflections  from  the  opposite  window.  For  these  mea¬ 
surements,  the  laser  beam  optics  were  designed  to  generate  a 
small  Gaussian  waist  at  the  focus.  A  low  laser  power  of  1 .0 
W  was  used  to  avoid  heating  of  the  dyed  fluid  in  the  very 
small  focal  volume  and  to  prevent  saturation.  The  plenum 
dye  concentration  was  kept  less  than  10  "  6  M.  Consequent¬ 
ly,  the  much  lower  concentrations  at  the  measuring  station 
did  not  significantly  attenuate  or  steer  the  beam. 

A  low  /  #  lens  was  used  to  collect  light  from  a  short 
segment  centered  at  the  waist  of  the  focused  laser  beam  (on 
the  jet  centerline,  at  x/d  =  100)  onto  a  photomultiplier 
tube,  effectively  yielding  single-point  concentration  values 
versus  time.  An  optical  low-pass  filter  eliminated  back¬ 
ground  laser  light,  transmitting  only  the  frequency-shifted 
fluorescence.  A  slit  spatial  filter  defined  the  length  of  the 
laser  line  segment  sampled.  The  slit  width  was  chosen  such 
that  the  sampling  volume  was  roughly  cylindrical,  about  80 
/zm  in  diameter  and  less  than  50/rm  in  length.  The  use  of  the 
slit  rather  than  a  pinhole  ensured  that  any  small  transverse 
beam  movements  did  not  alter  the  measurement  volume. 
The  long  Rayleigh  range  and  latitude  in  the  depth  of  field 


FIG  1  Experimental  apparatus. 
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minimized  the  effect  of  motion  in  either  of  the  two  horizontal 
directions.  Data  acquisition  was  computer-controlled.  The 
signal  amplifier  incorporated  a  three-pole  Butterworth  fil¬ 
ter,  with  a  cutoff  frequency  set  slightly  under  10  kHz.  The 
data  were  sampled  at  20  kHz  for  all  the  runs. 

Careful  consideration  has  been  given  to  the  reduction  of 
noise  in  these  measurements,  especially  since  we  wish  to  ex¬ 
amine  the  normalized  scalar  rms  fluctuations,  c'/  c.  As  de¬ 
scribed  previously, 17,1 8  we  utilized  Wiener  filter  tech¬ 
niques1920  in  processing  the  data.  The  issues  of  buoyancy 
and  attenuation  of  the  laser  beam  have  been  examined,  and 
are  discussed  in  Appendix  A.  We  conclude  that  neither  is  a 
factor  in  these  jets.  Finally,  an  important  consideration 
when  measuring  a  quantity  such  as  c’/  c  is  resolution,  since 
fluctuations  at  scales  smaller  than  the  measuring  volume  are 
low-pass  filtered.  It  is  therefore  necessary  to  determine  that 
either  the  smallest  scales  of  the  flow  are  resolved  or  that  the 
contribution  to  the  fluctuations  by  scales  smaller  than  the 
resolution  is  negligible.  We  discuss  resolution  concerns  later 
in  this  paper. 

III.  RESULTS 

For  our  discussions,  the  jet  Reynolds  number  is  defined 
as  Re  =  U„d /v,  where  Uit  is  the  jet  exit  velocity,  d  is  the 
diameter  of  the  jet  nozzle,  and  v  is  the  kinematic  viscosity  of 
the  working  fluids  (water). 

Our  measurements  indicate  that  the  scalar  field  of  a  high 
Schmidt  number  turbulent  jet  has,  in  some  sense,  not 
reached  an  asymptotic  state  for  Re  as  large  as  10  000  to 
20  000,  or  more.  This  is  apparent  in  the  behavior  of  the  scalar 
pdf ’s  with  varying  Re.  In  Fig.  2,  we  show  four  sample  pdf ’s, 
spanning  Reynolds  numbers  from  3000  to  24  000.  One 
hundred  and  fifty  histogram  bins  were  used  to  compile  each 
pdf.  Approximately  34  large-scale  structure  passages  were 
captured  at  the  lowest  Re,  and  over  125  at  the  highest. 

These  pdf ’s  display  a  decreasing  relative  concentration 
fluctuation  amplitude  and  imply  that  the  jet  is  becoming 
more  homogeneous  in  composition  with  increasing  Re.  It 
could  be  argued  that  this  may  be  attributable  to  worsening 
resolution,  as  Re  increases.  It  turns  out,  however,  that  the 


FIG  2.  Scalar  pdf's  for  varying  Re  (on  axis,  x/d  =  100) 


effect  is  not  due  to  resolution  difficulties,  as  will  be  addressed 
below. 

Another  manifestation  of  a  change  with  Re  is  in  the 
corresponding  scalar  power  spectra  (Fig.  3).  Here,  we  plot 
the  spectra  Ec  ( frs )  at  the  four  Reynolds  numbers,  normal¬ 
ized  by  the  estimated  large-scale  times, 

ta=6/UcI,  (1) 

and  the  local  mean  concentration  squared,  c2.  The  frequency 
axis  has  also  been  nondimensionalized  by  rb.  With  these  nor¬ 
malizations,  equal  ( equivalent )  length  scales  line  up  on  the 
fr6  axis,  regardless  of  Re.  The  large-scale  passage  frequen¬ 
cies,  as  independently  estimated  by  the  local  centerline  ve¬ 
locities  divided  by  the  width  of  the  jet,  are  seen  to  be  aligned 
at  /rA  =  1,  corresponding  to  the  length  scale  5.  So,  alterna¬ 
tively,  the  frb  axis  can  be  viewed  as  a  nondimensional,  in¬ 
verse  length  scale,  which  is  the  same  for  all  Re.  This  view  will 
be  useful  in  the  discussion  below. 

An  important  observation  is  that  these  spectra  do  not 
collapse  when  normalized  in  this  fashion.  The  central  ( mid¬ 
frequency)  portion  of  the  spectra  can  be  seen  to  bulge  up¬ 
ward  with  decreasing  Reynolds  number.  This  might,  per¬ 
haps,  be  understandable  given  the  low  Re  involved  and  the 
documented  reluctance  of  the  jet  spectrum  to  develop  a  con¬ 
stant  slope  region  until  at  much  higher  Re.  However,  in  ad¬ 
dition  to  the  shape  differences,  there  is  a  discrepancy  be¬ 
tween  the  low-frequency  levels  of  the  spectra.  If  the 
normalization  was  the  cause  of  these  discrepancies,  the  only 
“free"  parameter  would  be  rA,  the  estimate  of  the  large-scale 
time. 

Two  considerations  help  to  rule  out  attributing  the  var¬ 
iations  to  the  value  of  t„.  One  is  that,  since  t„  multiplies  the 
frequency  but  divides  the  spectrum,  changes  in  rA  shift  the 
normalized  spectra  along  lines  of  slope  —  1 .  This  does  not 
improve  the  collapse.  The  second  point  is  that  the  ratio  of  the 
constant  level  region  at  the  lowest  Re  to  that  at  the  highest 
Re  is  greater  than  2.  A  variation  in  rA  by  this  amount  is  not 
admissible  by  the  uncertainties  in  this  quantity.  Therefore, 
possible  uncertainties  in  normalization  cannot  explain  the 
variation  in  the  spectra  observed. 


log10(f  V 

FIG  3  Normalized  scalar  power  spectra  for  varying  Re  (on  axis. 

x/d  =  100) 


1158 


Phys  Fluids  A.  Vol  3,  No  5,  May  1931 


P  L.  Miller  and  P  E  Dimotakis 


1158 


Since  the  area  under  the  power  spectrum  is  equal  to 
c'2/  c2,  the  bulging  spectra  at  the  lower  Re  also  indicate  that 
c'/  c  decreases  with  increasing  Re.  It  is  significant  that  the 
area  of  greatest  contribution  to  this  effect  is  at  scales  unaf¬ 
fected  by  potential  resolution  difficulties,  namely,  from  the 
variation  at  the  lower  frequencies,  corresponding  to  time 
scales  of  the  order  of  the  large-scale  passage  time,  or  length 
scales  close  to  the  local  jet  diameter. 

We  calculate  the  normalized  scalar  rms  c'/  c  from  the 
data  and  plot  it  as  a  function  of  Re  ( Fig.  4 ) .  These  values  are, 
again,  from  data  collected  at  x/d  =  100,  with  the  exception 
of  one  point  (fifth  from  the  right),  which  is  from  data  re¬ 
corded  at  x/d  =  300.  The  resulting  c'/  c  can  be  seen  to  con¬ 
tinuously  decrease  with  increasing  Re.  This  is  consistent 
with  the  behavior  of  the  pdf ’s  and  the  spectra.  The  smooth 
line  in  Fig.  4  was  generated  by  fitting  the  measured  values  of 
c  /  c  using  a  function  of  the  form 

c‘/c  =  A+BRer.  (2) 

The  values  of  the  fit  parameters  obtained  were  A  =  0.20, 
B  =  550,  and  p  =  —  1 .02.  The  quality  of  the  fit  is  insensitive 
to  small  changes  (of  several  percent)  in  these  parameter  val¬ 
ues,  as  measured  by  the  mean  squared  error’s  response. 

We  note  that  the  nonzero  value  of  A  indicates  that  c'/  c 
may  be  approaching  an  asymptotic  value  of  about  0.20  at 
very  large  Reynolds  numbers.  We  will  return  to  this  point 
below. 

IV.  ISSUES  OF  RESOLUTION 

The  issue  of  resolution  for  these  measurements  is  crucial 
and  was  discussed  in  Miller  and  Dimotakis.17  To  reiterate, 
the  Kolmogorov’1  scale  can  be  estimated  directly  from  its 
definition  in  terms  of  the  mean  energy  dissipation  rate,  i.e., 

AK={v'/e)Wi.  (3) 

Using  the  result  of  Friehe  et  al.22  for  the  kinetic  energy  dissi¬ 
pation  rate  e  on  the  centerline  of  a  turbulent  jet,"  i.e., 


FIG.  4  Normalized  scalar  rms  fluctuations  for  varying  Re  (on  axis. 
x/d  =  100  and  300).  The  data  were  fit  with  a  function  of  the  form 
(  v  8  Re'\  The  parameter  values  determined  by  the  fit  are  A  =  0.20. 
B  =  550.  and  p  —  -  I  02 


e  =  48  ( Uo/d)[(x  —  x0  )/d  ]  ~  4,  (4) 

the  Kolmogorov  scales  for  these  measurements  are  found  to 
range  from  roughly  50  to  250  /im.  By  similarity  arguments, 
the  velocity  field  spatial  scale  where  the  action  of  viscosity 
will  become  important,  say,  A,.,  will  be  some  multiple  of  AK . 
Normalized  energy  spectra  are  found  to  break  from  a  con¬ 
stant  power  law  at  a  wave  number  k,.  such  that  k,AK  ~  1/8 
(e.g.,  Chapman23 ).  This  yields  an  estimate  of 

A,  =  Tr/k,.-~25AK.  (5) 

The  smallest  expected  scalar  diffusion  scale  A ,,  is  smaller 
yet  by  a  factor  of  Sc‘/2  (Batchelor24 ),  or  ~30  in  this  case, 
yielding  an  estimate  for  A,,  very  close  to  AK  (in  water). 
These  estimates  are  corroborated  by  the  gas  phase  experi¬ 
ments  of  Dowling  and  Dimotakis1 1  as  well  as  the  measure¬ 
ments  in  water  by  Buch  and  Dahm25  and  suggest  that  the 
smallest  diffusion  scales  may  have  been  resolved  in  these 
experiments,  at  least  at  the  lower  Reynolds  numbers. 

There  are  three  resolution  requirements  on  these  mea¬ 
surements:  spatial,  temporal,  and  signal-to-noise  (SNR).  If 
we  convert  our  spatial  resolution  estimate  to  the  units  of  Fig. 
3,  using  the  mean  centerline  velocities,  we  obtain  values  of 
about  fr„  =  1350.  Note  that  this  value  is  the  same  on  the 
nondimensionalized  frequency  axis  for  all  four  Re,  just  as 
with  the  large-scale  passage  times.  If  we  accept  that  the  high¬ 
est  resolvable  frequency  is  half  this  value,  we  estimate  our 
spectra  are  spatially  resolved  up  to  about  2.8  on  the  log  scale 
of  Fig.  3.  Thus,  at  the  lower  Reynolds  numbers,  the  spatial 
resolution  is  at  least  as  restrictive  as  SNR  resolution,  while 
for  the  highest  Re,  the  SNR  may  be  slightly  more  important 
in  determining  the  limit  of  resolution. 

To  estimate  the  impact  of  resolution  on  these  measure¬ 
ments,  we  utilize  a  common  classical  model  spectrum.  At 
sufficiently  high  Re,  one  might  expect  the  spectrum  to  have 
assumed  a  similarity  form.  We  approximate  this  as  a  con¬ 
stant  level  from  zero  frequency  to  the  large-scale  passage 
frequency,  a  region  of  constant  —  5/3  slope,  and  zero  be¬ 
yond  a  cutoff  at  f  t„  ( Fig.  5 ) .  The  constant  level  value  is 
defined  as 

K  =E.  ( /-,  </, r„ ) / c2rA .  (6) 


FIG  5.  Model  spectrum. 
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The  spectrum  is  normalized  here  to  the  variance  of  the 
scalar  fluctuations,  i.e., 


c'2 =2  r 

Jo 


Ec  ( /7 ) 


d(fr6). 


(7) 


which  relates  the  spectrum  directly  to  c'/  c.  We  integrate  the 
model  spectrum 


1  <2  rt»r„ 

K2  d(fT„ ) 

2  c2  Jo 

+  f,'V!(/r4)-w</(/r#) 


(8) 


to  obtain  the  result 

c'/c  =  tf{2  +  3[l  -  (/*//.  )1/,]},/J.  (9) 

On  the  jet  axis  and  at  fixed  x/d,  this  can  be  expressed  as  a 
function  of  Reynolds  number  by  substituting  for  /  and  / 
(cf.  Dowling 1,1 ),  yielding 

c'/c  =  *[2  +  3(1  -0.928  Re"  W1)]W2.  (10) 


Several  things  are  apparent  from  Eq.  ( 10).  First,  it  indi¬ 
cates  that  c' /  c  would  be  expected  to  increase  slightly  with 
Reynolds  number.  This  is  opposite  to  the  trend  exhibited  by 
our  data.  Also,  we  note  that  at  infinite  Reynolds  number  Eq. 
( 10)  predicts  a  finite  value  fore'/  cof  \/5  K.  To  evaluate  this 
expression  for  our  measurements,  we  fit  the  level,  low-fre¬ 
quency,  portion  of  our  highest  Re  spectrum  to  determine 
K \  yielding  *2~ 0.0078  and  an  asymptotic  value  of 
c'/c=:0.198.  This  value  of  K 2  agrees  well  with  the  corre¬ 
sponding  estimate  from  gas  phase  data, 1 1  while  the  asymp¬ 
totic  value  of  c'/c  agrees  with  the  suggested  infinite  Re  value 
of  0.20  obtained  by  fitting  the  data  [parameter^  in  Eq.  (2)]. 
Finally,  the  Reynolds  number  dependence  of  c'/  c  in  Eq. 
( 10)  is  extremely  weak.  For  Re  in  excess  of  10',  c'/c  has 
attained  over  99%  of  the  infinite  Re  asymptotic  estimate. 

The  model  spectrum  also  helps  demonstrate  the  effect 
resolution  limitations  would  have  on  the  resulting  measure 
'  f  c'/  c.  The  picture  is  much  the  same  as  in  Fig.  5,  but  we 
assume  that  the  measured  spectrum  is  resolution-limited  at 
some  nondimensional  frequency  frT„  (Fig.  6).  The  lost  con¬ 
tribution  toe'/  c,  represented  by  the  area  under  the  spectrum 
between /r,,  and /  r,,,  i.e., 

A  c'/c  =  c'/c  —  [c'/c],,  (11) 


is  equal  to 

Ac'/  c  =  \  3*  [  (/,//  )  ’  1  -  ( /,//  )-v>]  ,/-\  (12) 

while  the  ratio  of  the  resolution-limited  value  and  the  true 
value  is  given  by 

[cVc],  /2  +  3[1  -  f///  ):  '}  V  2 

c'/c  V 2  +  3 [  1  —  (/,//) 2/  ’]  / 

Here,  again,  the  ratio/,//  is  very  small  for  the  Reynolds 
numbers  of  interest.  As  a  result,  if  the  ratio/,//  is  also  kept 
small,  resolution  can  hardly  affect  the  measured  c'/c  at 
these  Re,  because  the  expression  in  Eq.  (13)  is  then  very 
close  to  unity.  Also  note  that  if  the  scalar  spectrum  is  ade¬ 
quately  described  by  the  model  spectrum  and  a  particular 
Reynolds  number  measurement  is  resolution  limited,  higher 


FIG.  6.  Resolution-limited  model  spectrum. 


Re  data  should  yield  the  some  (measured)  v->lue  of  [c'/  c], 
( not  less).  This  is  clear  because  the  measured  spectra  should 
then  have  the  same  shape  and  same  extent,  correspondingly 
yielding  the  same  values  of  [c'/  c]r.  If  a  Batchelor  A  1  re¬ 
gion24  was  incorporated  into  the  model  spectrum,  it  would 
not  alter  these  results  at  the  present  Schmidt  number  of 
about  10'  (see  Moffat21' ).  We  conclude  that  our  measure¬ 
ments  of  c'/  c  are  not  significantly  affected  by  resolution  diffi¬ 
culties. 

Finally,  recall  that  we  have  argued  that  we  may,  in  fact, 
be  resolving  the  smallest  diffusion  scales  in  the  flow,  at  least 
for  the  lower  Reynolds  number  flows  we  have  investigated. 
Specifically,  we  note  that  our  measured  spectra  do  not  ap¬ 
pear  to  support  the  k  1  Batchelor24  prediction  at  high  wave 
numbers  (frequencies).  This  observation  is  significant  at  the 
value  of  the  Schmidt  number  applicable  here,  since 
Scl/2~30.  This  would  have  yielded  a  k  1  portion  in  the 
spectrum  spanning  roughly  1.5  decades  in  frequency.  If  our 
assessment  of  the  resolution  of  our  measurements,  relative  to 
the  requisite  physical  scales,  is  correct,  then  the  resolution 
and  SNR  of  our  measurements,  at  least  at  the  lower  Reyn¬ 
olds  numbers,  are  sufficient  to  have  elucidated  this  behavior. 
As  regards  the  transport  of  scalars  characterized  by  low  dif- 
fusivity  (large  Sc)  in  turbulent  jets  in  the  Reynolds  number 
range  we  have  investigated,  we  are  forced  to  conclude  that 
they  do  not  display  a  k  1  Batchelor  spectrum. 

This  finding  should  be  examined  in  the  context  of  a  theo¬ 
retical  difficulty  with  the  A  1  spectrum.  In  particular,  the 
Batchelor  spectrum  does  not  appear  to  exhibit  the  correct 
limiting  behavior  at  high  Schmidt  numbers  ( holding  all  oth¬ 
er  flow  parameters  fixed).  As  Batchelor  had  noted,24  the 
variance  computed  using  the  A  1  spectrum  diverges  loga¬ 
rithmically  with  increasing  Sc.  This  divergence,  however,  is 
inconsistent  with  the  required  behavior  of  scalar  fluctu¬ 
ations,  which  must  remain  bounded,  irrespective  of  the  value 
of  the  Schmidt  number,  as  noted  recently.  Briefly,27  the  vari¬ 
ance  c’2  of  a  scalar  bounded  by  c  =  0  and  c  =  1  by  the  ini¬ 
tial/boundary  conditions  is  bounded,  in  turn,  by 

max{  c'2}  =  c(  1  —  c).  (14) 
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This  maximum  is  expected  in  the  limit  of  (identically)  zero 
diffusivity,  i/  —  0,  or  Sc  —  oc  if  all  other  flow  parameters  are 
held  fixed,  and  corresponds  to  a  limiting  pdf  of 

p(c)-(l  -c)S(c)  +c6(  1  -c), 
i.e.,  stirring  but  no  (molecular)  mixing  by  the  turbulent 
flow.  A  logarithmic  divergence  of  c'2  as  Sc  —  oo ,  as  estimat¬ 
ed  by  the  integral  over  the  k  1  scalar  spectrum,  is  inconsis¬ 
tent  with  the  bound  of  Eq.  ( 14). 


V.  CONCLUSIONS 

We  have  conducted  high-resolution,  single-point,  scalar 
measurements  on  the  centerline  in  the  far  field  of  an  axisym- 
metric,  turbulent  jet  at  high  Schmidt  number  Sc~  10',  over  a 
range  of  Reynolds  numbers  between  3000  and  24  000.  The 
jet  (Taylor)  microscale  Reynolds  number,  estimated  by" 

Re!  =  1.4  Re1/:,  (15) 

was,  therefore,  approximately  80  at  Re  =  3000  and  220  at 
Re  =  24  000.  Except,  perhaps,  for  the  lowest  values,  such 
microscale  Reynolds  numbers  are  regarded  as  high  by  some 
turbulence  criteria. :,i  Nevertheless,  our  experiments  indi¬ 
cate  that,  in  a  number  of  ways,  at  least  at  high  Schmidt 
numbers,  jet  mixing  has  not  attained  asymptotic  behavior  by 
Re=:2.4x  104.  In  particular,  pdf’s  and  spectra  of  c/c,  the 
normalized  jet  fluid  concentration,  are  found  to  depend  on 
Reynolds  number.  The  rms  concentration  fluctuations  c' /  c 
were  found  to  continuously  decrease  with  increasing  Re, 
over  the  Reynolds  number  range  examined. 

We  have  argued,  on  the  basis  of  our  spectral  measure¬ 
ments  and  the  use  of  a  simple  model  spectrum,  that  these 
results  are  not  an  artifact  of  inadequate  resolution.  In  partic¬ 
ular,  the  Reynolds  number  dependence  of  c'/  c,  as  estimated 
on  the  basis  of  the  model  spectrum,  should  have  exhibited 
the  opposite  behavior  to  that  observed.  The  limiting  values 
of  c  /  c  at  high  Re,  from  both  the  model  and  the  experimen¬ 
tal  estimates,  agreed  quite  well,  suggesting  that  c'/  c=:0.20 
may  represent  an  asymptotic  value  for  the  normalized  rms  of 
the  jet  fluid  concentration  fluctuations,  at  Sc  —  10',  on  the  jet 
axis.  Our  measurements  did  not  extend  to  jet  Reynolds 
numbers  high  enough  to  confirm  such  asymptotic  behavior 
directly. 

As  noted  earlier,  measurements  in  gas  phase  jets  at  com¬ 
parable  Reynolds  numbers' 1  suggest  that,  at  least  for  Sc  —  1, 
the  pdf  of  jet  fluid  concentration  does  not  exhibit  as  strong  a 
Reynolds  number  dependence  as  we  are  reporting  here,  if 
any  at  all.  This  difference  in  behavior  may,  of  course,  be  a 
manifestation  of  Schmidt  number  effects.  At  this  time,  we 
cannot  reconcile  the  disagreement  between  the  Reynolds 
number  dependence  of  the  present  results  and  the  liquid 
phase  “flame  length"  data  of  Weddell*  ( see  also  Hottel,  Ref. 
9).  and  Dahm  et  air"  In  particular,  on  the  basis  of  our  mea¬ 
sured  pdf’s  of  concentration,  we  would  conclude  that,  in 
liquid  phase  (Sc  -  10')  jets,  one  should  observe  decreasing 
flame  lengths  with  increasing  Reynolds  number,  at  least 
within  the  Reynolds  number  range  investigated  here  (see 
Appendix  B).  We  note,  however,  that  Kristmanson  and 
Danckwerts1"  found  the  flame  length  to  vary  with  Re  until 


about  twice  as  large  a  Reynolds  number  as  Weddell,  and 
Dahm  et  al. 

In  summary,  we  have  found  that  the  pdf ’s  of  the  jet  fluid 
concentration  become  narrower  with  increasing  Reynolds 
number.  A  measure  of  this  trend  is  provided  by  c'/  c,  the 
normalized  rms  level  of  scalar  fluctuations,  which  decreases 
as  the  Reynolds  number  increases.  These  measures  of  un- 
mixedness  suggest  a  more  homogeneous  turbulent,  liquid 
phase,  jet  fluid  concentration  field,  i.e.,  better  molecular 
mixing,  with  increasing  Reynolds  number.  Our  spectra  also 
evolve  with  increasing  Reynolds  number,  developing  a  con¬ 
stant  slope,  close  to  —  5/3,  at  the  higher  values. 

Finally,  we  did  not  observe  a  k  ~  1  spectral  behavior  at 
the  small  scales,24  even  though  the  predicted  range  of  valid¬ 
ity  of  this  regime  would  be  quite  large  at  this  Schmidt  num¬ 
ber  and  falls  within  the  resolution  range  of  our  measure¬ 
ments  at  the  lower  Reynolds  numbers  that  were 
investigated. 
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APPENDIX  A:  ISSUES  OF  BUOYANCY 

An  experimental  consideration  that  was  addressed,  in 
addition  to  resolution,  is  buoyancy.  There  are  two  cases  to  be 
considered:  buoyancy  from  a  jet  fluid/reservoir  fluid  density 
difference,  and  possible  heating  of  the  jet  at  the  measuring 
station  by  the  laser  beam.  Briefly,  a  nondimensionalized 
buoyancy  length  scale,'" 


( np>  v 

7  p.Ut,  V'-’ 

U  Pr ) 

\(P,  -P-T  )gd/ 

where  p ,  the  jet  fluid  density,  p ,  the  reservoir  fluid  density, 
(/„  is  the  jet  exit  velocity,  and  g  the  acceleration  of  gravity, 
may  be  derived  from  dimensional  analysis  (e  g.,  Dowling,"' 
Sec.  A.5).  This  parameter  may  then  be  compared  with  ex¬ 
perimentally  derived  values  to  determine  the  role  of  buoyan¬ 
cy  in  the  flow.  In  particular,  Papanicolaou  and  List14 "  find 
that  jets  are  momentum-dominated  if  x/Lh  is  less  than  one, 
while  Chen  and  Rodi"  cite  a  limit  of  roughly  0.5. 

We  will  estimate  the  maximum  anticipated  density  dif¬ 
ference  and  calculate  the  resulting  value  of  Lh/d.  The  den¬ 
sity  difference  may  arise  from  the  addition  of  the  heavy  dye 
to  the  jet  fluid,  or  from  a  temperature  mismatch  between  the 
jet  fluid  and  the  reservoir  fluid.  Let  us  consider  the  two  con¬ 
tributions  separately,  since,  it  turns  out,  they  are  of  different 
magnitudes.  The  fluorescein  dye  concentration  used  in  the 
experiments  was  roughly  3.3x10  7  M.  The  formula  weight 
of  disodium  fluorescein  is  376.28.  This  yields  an  estimate  for 
the  dye-induced  buoyancy  length  scale  [from  Eq.  (Al )  J  of 
over  19  600  jet  nozzle  diameters  at  the  lowest  Re.  For  each 
run,  the  jet  fluid  was  prepared  by  diluting  a  small  amount  of 
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concentrated  stock  solution  of  dye  with  sufficient  water,  tak¬ 
en  from  the  reservoir  shortly  before  the  run,  to  yield  the 
proper  jet  dye  concentration.  This  procedure  resulted  in  no 
observable  temperature  discrepancy  between  the  plenum 
fluid  and  the  reservoir  fluid.  Using  Eq.  ( A1 ),  the  stimated 
buoyancy  length  scale  resulting  from  a  temperature  differ¬ 
ence  of  0. 1  °C,  the  accuracy  of  the  thermometer,  is  1 500  jet 
nozzle  diameters.  While,  of  the  two  contributions,  the  possi¬ 
ble  influence  of  a  worst -case  temperature  mismatch  is  larger, 
the  resulting  buoyancy  length  is  15  times  the  distance 
between  the  nozzle  and  the  measurement  location  at 
x/d  —  100.  Therefore,  based  on  the  criterion  of  Chen  and 
Rodi,"  the  momentum-dominated  region  extends  for  a  dis¬ 
tance  that  is  over  seven  times  longer. 

A  second  way  buoyancy  may  arise  is  through  heating  of 
the  fluid  in  the  measuring  volume  via  absorption  of  laser 
power  by  the  dye.  The  approximate  residence  time  of  the 
fluid  in  the  measurement  volume  is  the  measurement  height 
(80^/m)  divided  by  the  calculated  mean  centerline  velocity 
at  the  lowest  Reynolds  number  (0.072  m/sec),  orO.OOl  sec. 
We  estimate  the  laser  attenuation  as  no  more  than  5%  across 
the  width  of  the  jet,  and  treat  the  jet  concentration  as  uni¬ 
form  in  the  radial  direction  for  the  purposes  here  ( not  too 
severe  an  approximation — cf.  Dahm12 ).  The  attenuation 
across  the  measuring  volume  (5x  10  7  cm')  would  then  be 
about  0.005%.  Assuming  a  worst  case  (all  power  attenu¬ 
ation  goes  into  heat),  the  resulting  temperature  rise  would 
then  be  approximately  0.02  °C.  While  the  approximation  of 
uniform  concentration  and  velocity  may  tend  to  underesti¬ 
mate  the  temperature  rise,  the  actual  laser  attenuation 
through  the  entire  jet  diameter  is  closer  to  3%,  or  so,  with 
less  than  half  of  the  attenuated  power  converted  into  heat. 
Therefore  this  estimate  is  conservative. 

A  second  point  concerning  the  heating  of  the  fluid  by  the 
laser  beam  is  that,  while  some  dilatation  of  the  fluid  occurs  as 
it  passes  through  the  laser  beam,  it  is  a  localized  effect.  Im¬ 
portantly,  there  is  no  global  buoyancy  force  to  affect  the 
large-scale  jet  structure  and  entrainment,  as  there  is  when 
the  jet  and  reservoir  densities  are  different.  Thus  it  is  not 
appropriate  to  utilize  criteria  such  as  those  of  Chen  and 
Rodi"  in  this  case,  since  they  apply  to  the  situation  with 
globally  different  densities.  The  primary  influence  here,  i.e., 
local  dilatation,  will  have  an  even  weaker  effect  on  the  flow 
for  a  given  temperature  rise.  The  dilatation  is  extremely 
small,  Ap/p~5x  10  ",  and  we  argue  that  its  influence  in 
our  flow  environment  is  also  insignificant. 

Having  addressed  the  issue  of  attenuation,  we  should 
also  consider  the  possibility  of  attenuation-caused  fluctu¬ 
ations  of  the  laser  beam  contributing  to  the  measured  values 
of  c'/  c  The  beam  attenuation  across  the  diameter  of  the  jet 
was  typically  3%,  or  less.  Fluctuations  arising  from  this  at¬ 
tenuation  may  be  completely  uncorrelated  with  the  scalar 
fluctuations  at  the  measurement  location,  in  which  case  the 
variances  of  the  two  contributions  add.  Alternatively,  the 
fluctuations  may  be  exactly  correlated,  in  which  case  the  two 
contributions  to  c7  c  add.  In  either  situation,  the  maximum 
influence  would  correspond  to  fluctuations  between  the 
minimum  (zero)  and  maximum  (  -  1.5%  )  attenuation  pos¬ 
sibilities.  or  a  maximum  correction  to  the  value  of  c7  c  of 


roughly  0.01 5.  This  would  result  in,  at  most,  a  7%  change  in 
the  measured  value.  This  should  be  compared  with  the  large 
percentage  decrease  in  the  measured  values  over  the  Reyn¬ 
olds  number  range  investigated  (Fig.  4).  Therefore,  even  in 
a  worst-case  scenario,  the  fluctuations  arising  from  the  at¬ 
tenuation  of  the  laser  beam  would  represent  only  a  small 
correction  to  the  magnitude  of  the  variation  we  have  mea¬ 
sured.  The  actual  situation  is  obviously  not  nearly  so  ad¬ 
verse,  and  we  conclude  that  the  attendant  contamination  of 
our  results,  as  a  result  of  beam  fluctuations,  is  extremely 
small. 

APPENDIX  B:  INTERPRETATION  OF  FLAME  LENGTHS 

Given  the  jet  fluid  concentration  pdf’s  p(c,x\ Re),  the 
“flame  length”  can  be  estimated  as  the  distance  required  for 
1  —  e  of  the  jet  fluid  (with,  say,  e  =  0.01 )  to  have  mixed  to 
the  proper  composition  (jet  fluid  fraction)  cr,  as  required  by 
the  chemical  reaction.  This  is  related  to  the  stoichiometric 
mixture  ratio  (equivalence  ratio)  d>,  by  the  expression 

c,=6,/(6,+  1),  (B1) 

where  6r  is  equal  to  the  parts  of  reservoir  fluid  that  must  be 
moleculariy  mixed  and  reacted  for  the  complete  consump¬ 
tion  of  one  part  of  (pure)  jet  fluid. 

Assuming  that  the  highest  values  of  the  jet  fluid  concen- 
trataion  occur  on  the  jet  axis,  to  estimate  the  flame  length 
corresponding  to  a  given  c,  and  e ,  we  need  two  additional 
assumptions.  First,  that  the  pdf  of  concentration  on  the  jet 
axis  is  a  function  of  c/  c  only,  i.e. 

p(c.x)  =/(c/c),  (B2) 

where  the  function  of f(c/  c )  may  vary  with  Reynolds  num¬ 
ber  (and  Schmidt  number),  but  is  taken  as  independent  of 
x/d.  Note  that  it  is  f(c/c)  that  was  plotted  under  this  im¬ 
plicit  assumption  in  Fig.  2.  This  is  the  usual  assumption, 
recently  corroborated  for  gas  phase  (Sc~l)  turbulent 
jets,11  with  supporting  evidence  for  liquid  phase  (Sc~  10') 
jets. 14  The  second  assumption  deals  with  the  dependence  of? 
on  x/d,  i.e.  (momentum-driven  turbulent  jets) 

c  =  xd  */(x  —  x„  ).  (B3) 

In  this  expression,  the  (pure)  jet  fluid  concentration  is  taken 
as  unity,  the  dimensionless  parameter  k  is  taken  as  a  constant 
( which  may.  however,  depend  on  Reynolds  number ) , 1 1  d  * 
is  the  jet  momentum  diameter,"  7 12  and*,,  is  a  virtual  origin 
(possibly  weakly  dependent  on  Reynolds  number  at  low 
Reynolds  numbers).  Combining  Eqs.  (B2)  and  (B3).  we 
have,  for  each  Reynolds  number 

p{c,x)  =/{[  (x  -  x„  )/Kd*]c}.  ( B4) 

For  a  given  e.  the  flame  length  L ,  may  then  be  estimated  as 
the  implicit  solution  of  the  equation 

J  p(c.L,)dc  =  e.  (B5) 

It  can  be  seen  that  the  computed  flame  length  is  determined 
by  the  behavior  of  the  pdf  at  high  values  of  c/  c.  Basically,  if 
/  (c/  ?)  is  the  pdf  of  c/c  at  a  Reynolds  number  Re;.  then  for 
c,/c>  1  (small  f)  and  f,  ( c,/c )  >/,  (c,/c),  we  will  have 
Li:  In  other  words,  and  as  perhaps  obvious  in  retro- 
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spect,  a  larger  c'/c  (less  well-mixed  jet)  implies  a  longer 
flame  length. 

The  inferred  flame  length  dependence  on  Reynolds 
number  made  above  assumes  that  any  deviations  of  the  x 
dependence  of p(cpc)  from  Eq.  (B2),  as  well  as  the  possible 
dependence  of  d  *,  k,  and  jc0  on  Reynolds  number,  are  not 
strong  enough  to  offset  the  relatively  large  dependence  of  L{ 
predicted  using  Eq.  (B4)  and  the  assumptions  outlined 
above.  We  note,  however,  that  the  validity  of  these  assump¬ 
tions  may  not  be  unassailable,  especially  at  the  high  Schmidt 
number  in  this  experiment.27 
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Abstract 

Results  from  experiments  performed  in  a  newly-constructed  High  Pressure 
Combustion  Facility  are  presented.  The  experiments  described  here  are  designed  to 
address  Reynolds  number  effects  on  turbulent  mixing  and  entrainment,  and  flame 
length  in  an  axisymmetric,  gas-phase  turbulent  reacting  jet,  in  particular,  over  a 
large  range  of  Reynolds  numbers. 


Introduction 

Flame  length  in  kinetically  fast  chemical  reactions  can  be  used  as  an  important 
measure  of  Reynolds  number  effects  and  can  be  related  to  the  rate  of  molecular 
mixing  in  turbulent  jets  (e.g.,  Broadwell  1982).  The  experiments  reported  here 
rely  on  a  new  technique  for  measuring  flame  length,  developed  for  the  purpose  of 
investigating  Reynolds  number,  as  well  as  other  effects  on  flame  length. 

The  facility  is  shown  schematically  in  Fig.  1.  The  main  apparatus  is  a  107  cm 
internal  diameter  by  183  cm  length  (1.64  m3)  pressure  vessel,  capable  of  operating 
at  pressures  from  0.1  atm  to  15  atm.  This  allows  an  independent  control  of  Reynolds 
number  over  2.2  decades  through  the  pressure,  while  holding  other  conditions  con¬ 
stant  (e.g.,  flow-time  scales,  nozzle  exit  conditions,  and  buoyancy  effects).  Buoyancy 
effects,  which  scale  as  Ap/poo,  where  Ap  is  the  average  change  in  density,  in  the  jet 
cone,  owing  to  heat  release  and  poo  is  the  ambient  density,  remain  constant  as  the 
pressure  is  varied  since  both  Ap  and  Poo  scale  directly  with  pressure. 
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FIGURE  1  Schematic  of  the  High  Pressure  Combustion  Facility. 

The  chemical  reactants  chosen  were  F2  and  NO,  diluted  in  N2,  for  several  rea¬ 
sons.  First,  these  reactants  are  hypergolic,  obviating  the  need  for  ignition  sources 
and  avoiding  flame  stability  issues.  The  resulting  chemical  kinetics  are  very  fast, 
ensuring  that  the  rate  of  chemical  product  formation  is  limited  by  the  turbulent 
mixing  process  and  not  by  the  chemical  kinetics  (high  Damkohler  number  regime), 
even  at  the  low  reactant  concentrations  employed  here.  Finally,  an  additional  at¬ 
tractive  characteristic  of  this  chemical  system  is  the  ability  to  produce  a  very  low 
adiabatic  flame  temperature  rise,  ATf  ~  7K  for  the  experiments  described  here, 
mitigating  buoyancy  effects  while  retaining  an  overall  reaction  rate  in  the  turbulent 
jet  that  is  kinetically  fast.  On  the  other  hand,  the  choice  of  these  reactants,  vis-a-vis 
safety  considerations,  dictates  a  strict  adherence  to  materials  properties,  as  well  as 
facility  assembly  and  maintenance,  and  experiment  operating  procedures. 

In  these  experiments,  the  jet  issues  from  the  bottom  of  the  pressure  vessel 
upwards,  through  a  2.5  mm  internal  exit  diameter  nozzle  that  protrudes  into  the 
vessel.  The  radial  profile  of  the  internal  nozzle  contraction  section  is  designed  as 
a  6th  order  polynomial,  optimized  to  avoid  Taylor-Gortler  vortices  in  the  concave 
parts  and  for  minimum  exit  boundary  layer  thickness.  The  outer  nozzle  profile  is  a 
paraboloid  of  revolution  in  an  attempt  to  cater  to  the  external  potential  flow  field 
that  is  induced  by  the  downstream  entrainment. 


The  jet  fluid  is  F2,  diluted  in  N2.  The  reservoir  fluid  is  NO,  also  diluted  in 
N2.  Partied  pressure  techniques  utilizing  intermediate  mixing  vessels  are  used  to 
accurately  produce  the  desired  jet  and  reservoir  concentrations.  The  jet  charge 
is  contained  in  a  smaller  volume  (0.043  m3)  reactant  tank  which  contains  a  teflon 
expulsion  diaphragm.  A  metered  sonic  N2  stream  is  fed  to  one  side  of  the  diaphragm, 
resulting  in  a  fixed  mass  flow  rate  of  jet  gas  through  the  nozzle.  The  reservoir 
pressure  is  held  constant  during  the  run  with  a  Camfiex  pressure  relief  valve.  See 
Fig.  1. 
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FIGURE  2  Schematic  of  jet  geometry  and  diagnostics  with  the  coordinate  axes  in¬ 
dicated. 


The  diagnostics  utilized  in  these  experiments  were  designed  to  give  an  accu¬ 
rate  measure  of  flame  length  using  constant  current  resistance  thermometers.  The 
thermometers  are  long  (58  cm  length),  thin  (25  pm  diameter),  Pt-10%Rh  wires, 
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stretched  across  the  jet  centerline  at  16  axial  locations,  from  xj  do  =  30  to  240. 
Here  x  is  the  streamwise  (jet-axis)  coordinate  and  do  is  the  nozzle  internal  exit  di¬ 
ameter.  The  wires  extend  well  beyond  the  edges  of  the  turbulent  jet  cone,  spanning 
the  heat  release  and  temperature  rise  regions,  at  all  downstream  locations.  The 
wire  streamwise  locations  were  selected  to  give  uniform  spacing  on  a  log10(x/do) 
axis.  This  spacing  was  suggested  in  view  of  a  conjecture  by  Dimotakis  that  the 
line-integrated  measurements  of  temperature  rise  in  the  flame  region  would  depend 
logarithmically  on  x/do  ■  Figure  2  shows  a  schematic  of  the  wires  with  the  coordinate 
axes  indicated. 

Each  wire  formed  one  of  the  legs  of  a  Wheatstone  bridge  circuit,  whose  output 
was  amplified  by  a  custom-built,  low-noise  signal  amplifier,  low-pass  filtered  by  a 
third-order  Butterworth  filter  with  a  knee  set  at  200  Hz,  sampled  at  500  Hz,  and 
digitized  to  12  bits.  The  wires  were  calibrated  by  allowing  the  interior  of  a  slightly 
overpressurized  vessel  to  equilibrate  in  temperature  overnight  and  allowing  the  pres¬ 
sure  to  decrease  over  a  few  seconds  in  the  morning.  The  resulting  resistance  change 
from  each  wire  was  recorded  during  the  adiabatic  cooling  process.  The  small,  com¬ 
mon  temperature  drop  experienced  by  each  wire  was  estimated  by  measuring  the 
internal  vessel  pressure  using  an  electronic  differential  pressure  transducer.  In  the 
course  of  a  chemically  reacting  run,  the  mean  temperature  rise,  determined  from 
the  change  in  total  resistance  of  each  wire,  was  estimated  in  this  fashion  to  within, 
approximately,  30  to  50  mK. 

The  reason  for  choosing  this  particular  measurement  technique  was  motivated 
by  the  scaling  laws  of  axisymmetric  free  turbulent  jets.  For  the  small  heat  release 
(temperature  rise)  employed  in  this  experiment,  the  local  wire  resistance  is  a  linear 
function  of  the  local  temperature.  As  a  consequence,  the  total  resistance  of  the 
cold  wires,  when  stretched  across  the  flame,  yields  an  accurate  estimate  of  the  line- 
integrated  heat  release,  i.e.,  temperature  rise,  resulting  from  the  chemical  reaction. 

At  the  end  of  the  flame  zone,  shown  schematically  in  Fig.  3,  the  F2  has  been 
completely  consumed.  No  heat  release  occurs  downstream  of  this  location.  Beyond 
the  flame  tip,  at  i  ~  If,  the  heat  released  by  the  reaction  upstream  is  simply 
diluted  by  the  subsequently  entrained  reservoir  fluid.  Therefore,  for  x  >  L(,  where 
Lf  is  the  flame  length,  temperature  excess  is  a  conserved  scalar. 


Assuming  that  the  flow  achieves  a  self-similar  temperature  profile  beyond  the 
flame  tip  (as  observed  by  Becker  &  Yamazaki  1978),  the  temperature  rise  for  x  >  L\ 
should  obey  the  jet  similarity  law  for  a  scalar  quantity  in  a  momentum-dominated 
jet.  This  similarity  law  can  be  derived  from  dimensional  arguments  (e.<?.,  Chen  & 
Rodi  1980),  has  been  verified  experimentally  ( t.g .,  Dowling  &  Dimotakis  1990),  and 
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FIGURE  3  Schematic  of  the  jet  flame  zone. 

where  A T(x,  y)  is  the  temperature  rise  at  (*,  y),  y  is  the  transverse  coordinate,  ATf 
is  the  adiabatic  flame  temperature  rise,  k  is  a  constant  (determined  by  experiment), 
d*  is  the  momentum  diameter  of  the  jet,  g(rj)  is  the  similarity  mean  profile  function 
(determined  by  experi  lent),  rj  is  the  similarity  variable  defined  as 
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with  xo  the  virtual  origin.  The  function  g(r})  has  a  maximum  value  of  1  at  rj  =  0 
and  goes  to  zero  outside  the  jet  edges  (rj  >  0.25).  The  momentum  diameter  d*  is 
considered  to  be  the  proper  normalization  for  the  downstream  coordinate  of  the  jet 
(e.y.,  Thring  &  Newby  1953,  Avery  &  Faeth  1974,  and  Dahm  &  Dimotakis  1987), 
and  is  defined  as 


dT  = 


2  mp 
y/npooJo 


(2) 


where  rho  and  Jq  are  the  nozzle  exit  mass  and  momentum  fluxes,  find  pe «  is  the 
reservoir  gas  density. 
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Taking  the  line  integral  of  Eq.  la  along  y,  at  constant  x  (along  a  wire),  and 
realizing  that  g(i 7)  goes  to  zero  before  y  =  ±  Lm/2,  yields 


_1_  [L-n  A T(x,y) 
iw  J-Lw/2  AT, 


dy 


■  *  (£)  /Ij(’,)d’' 

^  fa(»)  ,  (3) 

where  Lw  is  the  span  (length)  of  the  wire. 


As  noted  above,  the  measurement  performed  by  each  wire  is  the  line  integral 
of  the  temperature  excess  at  constant  x.  This  is  given  by 


1  /-W2 

7-  /  AT(x,y)dy  . 
J~Lml 2 


(4) 


This  quantity  is  analogous  to  the  product  thickness  used  for  the  shear  layer  (e.g., 
Dimotakis  1989).  In  particular,  the  product  thickness  6p(x),  normalized  by  Lw, 
will  be  deiSned  as 


^p(x) 

Iw 


~  iw  y_oo  at, 


dy  . 


(5) 


The  product  thickness  6p(x)  is  equal  to  the  width  of  a  top-hat  temperature  profile  of 
height  AT,  that  yields  the  same  area  as  that  found  under  the  A T(x,  y)  temperature 
rise  profile  for  that  x  location.  Equations  3  and  5  are  equivalent  to  Eq.  4  divided 
by  AT,. 


The  time-averaged,  line-integrated  mean  temperature  rise  measured  at  each 
wire  beyond  the  flame  tip,  *.e.,  for  x  >  Lf,  is  expected  to  asymptote  to  a  constant 
value.  Accordingly,  the  fine  integral  will  rise  to  that  asymptotic  value  as  dictated 
by  the  cumulative  mixing  and  chemical  product  formation  up  to  that  location  for 
values  of  x  <  Lf. 


Results  and  Conclusions 


Figure  4  shows  a  sample  of  the  raw  data  from  the  16  wires  spanning  x/dQ  =  30 
to  240,  along  with  two  reference  wires  placed  outside  of  the  jet  cone  at  x/do  =  0 
and  240,  and  the  scaled  velocity  trace.  Starting  from  the  bottom  trace  and  working 
upwards,  the  lowest  trace  is  the  reference  wire  at  x/do  —  0  (partially  obscured  by 
the  x-axis  of  the  figure),  then  the  16  measurements  from  x/d0  =  30  to  240,  the 
other  reference  wire  at  x/do  =  240,  and,  finally,  the  scaled  velocity  trace.  The 
reference  wire  at  x/do  =  0  measures  the  temperature  of  the  ambient  reservoir  fluid 
and  the  reference  wire  trace  at  x/d0  =  240  indicates  the  arrival  of  hot  recirculated 
products  at  the  last  measuring  station,  thus  determining  the  temporal  extent  of 
uncontaminated  data.  The  scaled  velocity  trace  indicates  the  time  required  to 
establish  steady  state  flow. 


Each  channel  offset  by  0.75 K 


FIGURE  4  Sample  traces  of  cold  wire  raw  data  for  <f>  =  15  and  A  7}  =  7K.  Vertical 
dashed  lines  mark  the  averaging  interval  used  for  the  corresponding  data 
in  Fig.  5.  See  text  for  details. 
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In  the  typical  data  in  Fig.  4,  note  the  small  delay  between  the  time  the  jet 
velocity  has  reached  its  steady  value  and  the  time  the  temperature  rise  is  registered 
in  the  first  wire  at  x/do  =  30.  This  time  corresponds  to  the  time  required  to  flush 
the  jet  plenum  of  its  original  contents  of  reservoir  fluid  and  charge  it  with  the  F2- 
bearing  jet  fluid.  Note  also  the  increasing  delay  required  for  the  reaction  head  to 
reach  the  wires  at  increasing  x/do.  The  first  vertical  dashed  line,  at  ti  ~  3.3 sec, 
marks  the  beginning  of  the  time  interval  over  which  the  data  from  each  wire  are 
to  be  averaged.  The  second  vertical  dashed  line,  at  t2  c*  13.6  sec,  marks  the  end  of 
the  averaging  interval,  prior  to  the  arrival  of  the  recirculating  hot  product  gases,  as 
registered  on  the  top  reference  wire  at  x/do  =  240,  outside  the  turbulent  jet  cone 
(second  trace  from  the  top).  Note  that  the  absolute  span-averaged  temperature 
change  registered  by  the  first  wire  is  approximately  0.5  K.  Small  changes  in  the 
vessel  temperature  in  the  course  of  the  run,  as  registered  by  the  reference  wire  at 
x/do  =  0,  are  subtracted  from  the  temperature  registered  by  each  wire  and  the 
difference  for  each  wire  is  averaged  over  the  selected  (<1 ,  #2 )  time  interval  for  each 
run.  As  can  be  seen  from  the  time-trace  of  the  wire  at  x/do  —  30,  the  time-averaged 
mean  temperature  rise  is  determinable  to  a  very  good  fractional  accuracy. 


The  run  recorded  in  Fig.  4  was  for  a  stoichiometric  mixture  ratio  of  <f>  =  15 
and  an  adiabatic  flame  temperature  rise  of  ATf  =  7  K.  The  stoichiometric  mixture 
ratio  4>  represents  the  mass  of  reservoir  fluid  that  must  be  entrained,  mixed  on  a 
molecular  scale,  and  reacted  to  completely  consume  a  unit  mass  of  jet  fluid.  It  is 
defined  by 

/  _  (1W>N0) 

”  (YrJYvo)„  ’ 

where,  for  a  mixture  of  N  species,  the  mass  fraction,  Yi,  is  defined  as 


y,  2 


rrij 

N 

E  m« 

1=1 


(7) 


with  m,  equal  to  the  mass  of  the  ith  species  in  the  mixture.  The  resulting  averaged 
data  points  from  the  run  in  Fig.  4  are  depicted  in  the  corresponding  <j>  =  15  trace 
in  Fig.  5. 


Data  for  a  number  of  stoichiometric  mixture  ratios  4>  and  an  adiabatic  flame 
temperature  rise  of  ATf  =  7  K  were  processed,  as  described  above,  to  produce 
the  plot  depicted  in  Fig.  5.  The  data  correspond  to  the  nondimensional  product 
thickness  (Eq.  5)  and  sire  plotted  versus  log10(x/<f*)  for  the  four  stoichiometric 
mixture  ratios  4>  indicated.  The  solid  line  for  each  set  of  data  is  a  compound  curve 
fit  with  a  linear  least  squares  fit  to  the  ramp  region  and  a  cubic  spline  fit  to  the 
knee  region,  with  the  appropriate  matching  of  value,  slope,  and  curvature  at  the 
crossover.  The  flame  length  Lf  will  be  defined  as  the  99%  of  maximum  6p/Lw  point, 
in  the  same  spirit  as  the  definition  of  a  99%  boundary  layer  thickness.  The  expected 


1.4  1.6  1.8  2.0  2.2  2.4  2.6 

log10(x/cJ*) 

FIGURE  5  Sample  plot  of  product  thickness  vs.  log10(x/d*)  for  several  <j>' s  and  an 
adiabatic  flame  temperature  rise  of  AT{  =  7  K. 

increase  in  the  integrated  temperature  rise  is  observed,  as  well  as  the  asymptotic 
value  at  the  end  of  the  flame  tip.  The  constant  temperature  level  me  isur  1  by  the 
wires  beyond  the  flame  tip  is  proportional  to  <j>,  as  can  be  shown  from  simple  scaling 
arguments.  The  results  display  a  linear  behavior  in  the  flame  region,  when  plotted 
in  the  semi-logarithm i c  coordinates  suggested  by  the  choice  of  wire  spacing,  ai.i 
corroborate  the  conjecture  by  Dimotakis  mentioned  above. 

The  results  of  additional  experiments  will  be  reported  in  the  near  future  that 
address  the  flame  length  dependence  on  Reynolds  number,  the  transition  of  the  jet 
from  a  momentum-dominated  to  a  buoyancy-dominated  regime,  and  the  effects  of 
finite  rate  kinetics  (Dahmkohler  number)  on  mixing  in  the  turbulent  jet. 
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can  depend  on  the  elements  used  to  cover  the  set  (see.  for  ex¬ 
ample.  Feder  [1988.  p.  14] ).  In  the  Mandelbrot  1977  and  1982 
discussions,  the  two  quantities  arc  used  interchangeably,  by  con¬ 
sidering  cases,  as  he  says,  for  which  this  equivalence  can  be 
accepted 

As  Mandelbrot  notes  |I982.  pp.  357-360],  the  basic  mathe¬ 
matical  ideas  can  be  traced  to  earlier  proposals  and  discussions. 
Equation  (2)  was  essentially  used  by  Ponttjagin  and  Schnirelman 
in  the  30s.  and  the  quantity  D  defined  by  Eqs.  ( 1 1  and  (2)  can  be 
identified  with  the  Kolmogorov  capacity  (Kolmogorov  and 
Tihomirov.  1959).  or.  as  it  was  subsequently  dubbed,  the  box¬ 
counting  dimension  There  can  be  no  doubt,  however,  that  it  was 
Mandelbrot  who  realized  the  importance  of  these  ideas,  brought 
them  to  the  attention  of  a  community  outside  that  of  pure  math¬ 
ematics,  and  argued  that  they  should  be  expected  to  prove  appli¬ 
cable  and  useful  in  the  description  of  many  physical  phenomena. 
A  considerable  body  of  research  has  appeared  on  the  record  in 
response  to  his  proposals  relating  to  a  spectrum  of  phenomena  too 
numerous  to  trace. 

In  the  context  of  many  of  the  phenomena  proposed  by  Man¬ 
delbrot  to  be  potentially  dcscribable  by  the  fractal  power  law 
coverage  similarity,  however,  one  must  require  that  Eqs.  ( 1 )  and 
(2)  hold  for  X  not  only  strictly  positive,  but  over  some  non  van¬ 
ishing  range  of  scales  A.  i.e..  for 

0  <  Xmm  s  A  S  Xmj,.  (4) 

This  is  a  necessary  assumption  if  one  is  to  derive  many  of  the 
benefits  of  this  analysis.  It  will  also  prove  important  on  another 
score . 

Specifically,  to  return  to  the  basic  proposal,  it  is  important  to 
note  that  as  regards  the  description  of  physical  entities  in  which 
X  has  units.  Eq.  ( I )  is  not  dimensionally  correct ;  only  a  dimen¬ 
sionless  group  can  be  equal  to  a  pure  number.  Of  course,  one 
could  argue  that  a  proportionality  constant  C  is  implicit  in  Eq. 
(1).  i.e.. 

MX)  =  CX  "  (5) 

Such  a  constant  would  be  dimensional,  however,  with  units  as 
dictated  by  those  of  X  and  the  value  of  D.  i.e..  C  =  fn<Dl  In 
using  Eq  (5)  as  an  implicit  definition  of  D,  in  other  words,  in 
which  D  is  the  independent  variable,  the  "constant''  of  propor¬ 


tionality  C  is  not  a  constant  at  all  and  hides  important  scaling 
information. 

To  illustrate  the  point,  consider  a  set  representing  a  surface  in 
a  d-dimensional  (Euclidean)  space,  for  example,  fixed  threshold 
crossing  of  a  scalar  time  trace,  or  linear  image  (d  =  1 );  isoscalar 
contours  of  a  scalar  field,  e.g..  d.v.v)  =  clh.  in  a  planar  image  id 
=  2);  or  isoscalar  surfaces,  e.g..  cix.y.:)  =  c,h.  of  a  scalar  field 
in  three  dimensions  (d  =  3).  Consider  a  volume  given  by  L‘‘.  in 
each  case,  where  L  is  the  (linear)  extent  of  the  space  spanned  by 
the  data.  The  total  number  of  X-sized  elements  that  can  fit  in  this 
space  is  given  by  (L/AV'.T  The  total  coverage  Nj  X )  of  the  set  can 
then  be  expressed  as  the  fraction  of  the  number  of  X-elements 
covering  the  set.  times  the  total  number  of  X-clements  in  the 
^-dimensional  (data)  space,  i.e.. 

Nj  (X)  =  Fd  (X)  ^ V  (6) 

The  proposal  that  Nd(\)  *  \~Dj  is  then  equivalent  to  a  depen¬ 
dence  of  the  coverage  fraction  on  X.  i.e..  the  probability  that  a 
A-element  covers  a  portion  of  the  set.  given  by 

Fd(\)  *  (7) 

Since  Fd  (X)  is  a  probability  and  therefore  dimensionless.  Eq.  7 
requires  that  X  must  be  scaled  by  some  like-dimensioned  quan¬ 
tity.  or  characteristic  scale.  A0.  i.e., 

/  X  \J~Dj 

fj(A)  =  Fj0  (8) 

We  can  combine  this  with  Eq.  6  to  yield  what  is  implicit  in  the 
original  equation,  namely. 

(j)' 

In  this  expression,  the  now  dimensionless  constant  F^,,  is  equal  to 
the  required  fraction  of  set-covering  elements  of  size  X  =  A„. 
provided,  of  course,  the  characteristic  scale  A„  falls  within  the 
(Xmin.Xm„)  interval  in  which  Eq.  9  applies.  While,  in  principle, 
the  factor  F M  could  be  absorbed  in  a  redefined  X„.  it  is  useful  to 
retain  it;  the  characteristic  scale  A„  might  well  be  definable  on 
some  other  physical  basis  quite  independently  of  a  fit  to  the 
coverage  of  the  set  of  interest  using  Eq  9. 


t  Barring  edge  effects,  which,  however,  can  be  eliminated  in  any  of  a  number  of 
wavs — for  example,  bv  considenng  A  s  that  are  exact  binarv  submuluplcs  of 

L.  i.e..  X.  «  LT. 
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As  an  aside,  it  might  appear  that  Eq.  (3).  the  Hausdorff- 
Besicovitch  expression  for  D.  is  not  subject  to  the  same  diffi¬ 
culty.  Strictly  speaking,  it  is  dimensionally  correct  as  it  stands. 
Nevertheless,  one  would  prefer  not  to  be  using  a  measure  that  is 
dimensional  (with  units  that  depend  on  the  value  of  D).  or  to  be 
taking  the  limit  with  a  dimensional  A  — *  0.  to  determine  the 
critical  value  of  D.  Instead,  one  would  prefer  a  measure  based  on 
a  scaled  k.  e  g..  .<t(A/A0;  D).  with  D  computed  as  the  critical 
exponent  in  the  limit  of  A/A0  — »  0.  Again,  in  other  words,  there 
is  a  need  for  a  characteristic  scale  A„.  even  if  not  as  manifestly  as 
with  the  box-counting  coverage  expression  of  Eq.  ( I ). 

These  observations  have  an  interesting  consequence. 

Returning  to  Eqs.  (1),  (5).  or  (9).  we  see  that  they  can  be 
expected  only  to  apply  to  processes  that  observe  a  characteristic 
scale  A0.  The  resulting  coverage  function  MX)  could  then  depend 
strongly  on  this  characteristic  scale,  or  not.  depending  on  the 
value  of  the  dimension  exponent  D.  More  importantly,  however, 
the  coverage  of  sets  generated  by  processes  that  do  not  possess 
such  a  characteristic  scale  cannot  be  expected  to  be  described  by 
these  equations. 

Turning  to  the  more  specific  question,  an  area  in  which  there 
have  been  many  attempts  to  apply  power  law  fractals  is  ".  .  .the 
grand  chapter  of  Physics,  the  study  of  turbulence  in  fluids" 
(Mandelbrot.  1977,  p.  1 45 ( — the  dynamics,  in  other  words,  of 
high  Reynolds  number  fluid  motion.  Turbulent-flow -generated 
sets,  especially  at  high  Reynolds  numbers,  are  known  to  exhibit 
a  host  of  similari'y  properties.  These  can  be  expected  to  hold  over 
a  large  range  of  space  and  time  scales  that  are  bracketed  by 

~  S.  (10a) 

the  outer  flow  scale,  and 


the  inner  Kolmogorov  1 1 94 1  a |  scale,  in  which  v  is  the  kinematic 
viscosity  and  i  is  the  mean  dissipation  rate  of  the  kinetic  energy 
per  unit  mass.  At  high  Reynolds  numbers  Re.  the  bounds  of  this 
(inertial)  range  arc  in  the  ratio 

7^~p~Rc,J.  (ID 

This  potentially  very  large  similarity  range  would  seem  to 
provide  a  promising  arena  for  the  Mandelbrot  fractal  proposals. 
Sets  generated  by  turbulent  flow  processes  that  have  been  ana¬ 
lyzed  in  this  light  include  selected  velocity  component  time 
traces,  estimates  of  the  local  energy  dissipation  and  its  fluctua¬ 
tions,  temporal  as  well  as  spatial  scalar  fluctuations  and  isoscalar 
surfaces,  etc.  The  recent  review  by  Srccmvasan  |199l|  can  be 
consulted  for  a  discussion  and  references. 

Let  us  examine  this  proposal  in  view  of  some  of  the  similarity 
properties  in  the  inertial  range  of  scales  and  the  preceding  dis¬ 
cussion. 

In  his  first  and  second  similarity  hypotheses.  Kolmogorov 
( 1 94 1  a |  proposed  that  the  only  dynamically  important  quantity  at 
a  scale  X  =  2tr Ik  within  the  range  of  Eq.  (4).  with  the  identifi¬ 
cations  ip  Eqs  ( 10).  is  the  mean  dissipation  rate  i.  This  proposal 
immediately  leads  to  the  celebrated  Kolmogorov  -  V<  spectrum 
by  dimensional  analysis,  i.e.. 

£<fc)  =  CKe:  V-  '.  (12a) 

where  CK  is  a  (dimensionless)  constant  of  proportionality  of  order 
unity  The  right-hand  side  is  the  only  group  that  can  be  formed 
with  «  and  k  that  has  the  requisite  dimensions. 

Physically,  the  first  and  second  Kolmogorov  hypotheses  are 
equivalent  to  assuming  that  the  dynamics  of  the  turbulent  strain¬ 


ing  motions  that  take  an  eddy  of  a  scale  X,  to  a  scale  A:.  for 
&  ^  A,  >  Xt  ^  Xk. 

depend  only  on  the  local  scale  and  therefore  only  on  the  ratio 
A, /A,.  In  particular,  they  are  equivalent  to  assuming  that  the 
dynamics  do  not  depend  on  the  ratio  of  the  local  scale  X  to  the 
outer  scale  5.  the  inner  scale  XK  or  any  other  scale. 

In  subsequent  refinements  by  Kolmogorov  ( 1962],  Oboukhov 
1 1962].  and  others,  the  kinetic  energy  dissipation  rate  e  =  «(x,r) 
was  accommodated  as  the  local,  intermittent  field  that  it  is.  This 
introduced  the  outer  scale  8  as  a  scale  weakly  admissible  in  the 
dynamics,  resulting  in  small  corrections  to  the  -Vs  spectrum. 
Specifically,  with  a  log-normally  distributed  e(x,r).  the  energy 
spectrum  becomes 

Elk)  =  CK  €-  '  (kh)-'x'>  k-*\  (12b) 

where  p.  is  the  intermittency  exponent,  estimated  to  be  in  the 
range  of  0.2  <  p  <  0.5  (e.g..  Monin  and  Yaglom  1 1 975 .  p. 
642]).  Significantly,  the  difference  between  Eqs.  ( 12a)  and  ( 12b) 
is  too  small  to  be  discernible  experimentally.  At  least  by  this 
measure,  we  may  conclude  that  the  Kolmogorov  [  1941a]  hypoth¬ 
esis  of  scale  independence,  that  leads  to  the  -  Y>  spectrum  law  , 
must  very  nearly  be  right. 

This  conclusion  has  consequences,  in  light  of  the  preceding 
discussion  of  the  need  for  a  characteristic  scale  to  exist,  if  the 
power  law  fractal  similarity  is  to  describe  the  coverage.  In  the 
case  of  sets  with  scales  in  the  inertial  range  and  the  (near)  absence 
of  a  characteristic  scale  X„.  we  cannot  expect  Eq.  (9)  to  apply  if 
the  fractal  dimension  exponent  Dy  for  coverage  of  three- 
dimensional  sets  is  to  be  assigned  values  in  the  range  of  /),  = 
2.36  t  0.05.  as  was  suggested,  for  example,  by  Sreenivasan 
(1991.  p.  556].  This  is  also  a  value  close  to  £>,  =  Vs  that  Man¬ 
delbrot  suggests  "is  reducible  to  the  Kolmogorov  theory  and  the 
empirical  spectra. "+  which  he  also  argues  corresponds  to  a  value 
of  D2  =  Vs  for  planar  sections  (Mandelbrot  1977.  p.  52]. 

We  note  here  that,  implicit  in  the  assignment  of  both  the  re¬ 
ported  experimental  estimate  of  D,  =  2.36  ±  0.05  of 
Sreenivasan.  as  well  as  the  correspondence  between  the  value  of 
Dy  =  Vs  and  D:  =  V<  proposed  by  Mandelbrot  as  applicable  to 
isoscalar  surfaces,  is  the  relation 

d |  -  Dd<  =  d2  -  Dd_.  (13) 

where  d,  and  d:  denote  the  topological  dimension,  e  g  .  D,  =  Dz 
+  I  and  Dy  =  Dt  +2.  as  was  asserted  by  Sreenivasan  and 
Mencveau  (1986.  Sec.  2.2] — recall  also  Eqs.  (8)  and  (9).  Spe¬ 
cifically.  in  the  case  of  the  experimental  values  reported,  the 
value  of  Dy  cited  was  inferred  from  measurements  of  one¬ 
dimensional  (d  =  I )  or  two-dimensional  (d  =  2)  data,  assuming 
Eq.  (13).  A  suggestion  to  that  effect  had  previously  been  made  by- 
Mandelbrot.  whom  Sreenivasan  and  Mencveau  (1986)  cite  as  a 
reference  for  this,  and  who  says  that  "it  is  almost  sure-  "  that 
the  intersection  of  a  three-dimensional  surface  with  a  plane  yields 
a  dimension  Dz  =  Dy  -  I  (Mandelbrot.  1982.  p.  135],  espe¬ 
cially  if  the  three-dimensional  surface  of  interest  is  irregular 
(ibid.,  p.  136).  The  latter  proviso  would  of  course  be  satisfied  in 
the  case  for  an  isoscalar  surface  in  turbulence,  for  example  (see 
also  the  discussion  in  Mandelbrot  (1976.  p.  124;  1977.  p.  52 1 
cited  above).  We  will  return  to  this  point  below. 

Accepting  Eq.  (13)  at  face  value  for  now  .  we  can  argue  that 
such  values  for  D,  for  an  isoscalar  surface,  for  example,  are  too 
far  from  the  topological  dimension  of  a  three-dimensional  object, 
i.e..  D,  =  3.  and.  despite  Mandelbrot's  assertion,  that  the  re¬ 
sulting  dependence  of  the  coverage  function  /V,(A)  on  whatever 

th  is  not  clear  whether  Mandelbrot,  in  his  reference  here  to  ‘the  Kolmogorov 
theory  and  the  empirical  spectra.'*  is  thinking  of  the  original  <  1941a)  argu¬ 
ments  leading  to  the  -  Vi  spectrum  |Eq  ( 12a).  or  the  revised  ( 1962)  proposals 
that  lead  to  Eq  (12b)| 


Nonlinear  8cl«nct  Today 


1991  Vol.  1,  No.  2 


29 


characteristic  scale  Xo  is  used  in  Eq.  (9)  is  too  strong  to  be 
consistent  with  turbulent  flow  inertial  range  similarity.  This  con¬ 
clusion  can  be  reached  by  considering  the  properties  of  an  iso¬ 
scalar  surface  and  its  proposed  (power  law)  fractal  coverage 
N3(X),  assuming  that 

Xo*8,  (14) 

as  one  might  in  view  of  the  turbulent  cascade  process  from  larger 
to  smaller  scales.  If  we  consider,  for  example,  two  gedanken 
experiments  in  the  same  type  of  flow,  arranged  to  generate  the 
same  inner  scales  XK,  but  with  outer  scales  5,  and  h2  differing  by, 
say,  a  factor  of  10.  and  compare  the  coverage  N,(X)  at  a  fixed 
scale  X,  say,  for  X  =  10XK  «  The  analysis  of  such  a 

comparison  is  beyond  the  purposes  of  the  present  brief,  informal 
discussion  and  is  not  included  here. 

This  argument  was  recently  noted  in  Miller  and  Dimotakis 
[1991],  who  did  not  find  Eq.  (1)  [or.  equivalently.  Eq.  (9)] 
applicable  to  measurements  of  temporal  isoscalar  intersections  at 
a  fixed  point,  or  to  spatial  measurements  derived  from  linear 
images,  in  the  far  field  of  turbulent  jets  where  it  had  previously 
been  reported  to  apply  (e.g..  Sreenivasan  and  Meneveau  [  1986]). 
In  a  note  added  in  proof,  Sreenivasan  [1991]  suggested  that  the 
discrepancy  could  be  attributable  to  the  fact  that  Miller  and  Di¬ 
motakis  analyzed  temporal  data  As  noted  above,  however,  not 
only  had  such  behavior  been  reported  previously  for  temporal 
data,*1  but  in  their  discussion  Miller  and  Dimotakis  included 
measurements  and  the  analysis  of  temporal,  spatial,  as  well  as 
spatio-temporal  (streak  image)  data,  with  similar  results. 

Briefly,  to  cover  a  total  of  A/T  transitions  in  a  one-dimensional 
record  of  length  L.  we  must  find 


for  X  — >  0 
for  X  >  8. 


(15) 


The  small-scale  limiting  behavior  can  be  understood  by  appreci¬ 
ating  that,  once  the  covering  element  scale  X  has  decreased  below 
the  shortest  separation  between  threshold  transitions  in  the  scalar 
trace,  the  total  number  of  elements  of  size  X  required  to  cover 
those  transitions  becomes  equal  to  the  total  number  NT  of  tran¬ 
sitions  in  the  record  of  length  L 

The  large-scale  limiting  behavior  can  be  understood  by  appre¬ 
ciating  that  every  covering  element  of  size  X  much  larger  than  the 
outer  scale  8  will  cover  transitions  and.  therefore,  the  total  num¬ 
ber  of  elements  required  to  cover  those  transitions  must  tend  to 
the  total  number  of  elements  Z./X  that  can  fit  in  the  record.  In 
other  words,  we  must  have  F ,<X)  — *1  for  X  8  (cf.  Eq.  6).  It  is 
also  easy  to  appreciate  that,  for  typical  traces  derived  from  tur¬ 
bulent  flow,  the  number  of  elements  A/,(X)  needed  to  cover  the 
corresponding  sets  would  only  gradually  approach  the  limiting 
values  in  Eq  (15).  for  small  and  large  X. 

In  view  of  the  above  and  Eq.  (2).  we  must  have,  for  one¬ 
dimensional  isoscalar  data,  that 


fO  for  X  <  XB 
1 1  for  X  >  8. 


(16) 


where  XB  is  the  Batchelor  |I959|  scalar  diffusion  scale  (below 
which  the  isoscalar  surface  can  be  regarded  as  differentiable),  and 
that,  therefore,  the  exponent  D,  cannot  be  treated  as  a  constant. 

While  it  may  be  a  minor  point,  it  should  be  noted  that  if  the 
measurements  are  alias-free  with  some  margin,  i.e..  if  they  have 
been  sampled  at  a  frequency  comfortably  higher  than  the  Nyquist 


Figure  1:  Coverage  W,(X)  derived  from  temporal  data  of 
jet  fluid  concentration  crossings  of  the  local  c, 
measured  on  the  axis  of  a  turbulent  jet  at  xid  =  100 
and  a  Reynolds  number  of  1.15  x  10<.  The  scale  X  here 
is  converted  to  meters.  Data  from  Miller  and  Dimotakis 
[1991]. 


frequency,  and  sufficiently  noise-free.  MX)  should  tend  to  a  con¬ 
stant  at  the  smallest  values  of  X  resolved  by  the  data  (for  the 
reasons  stated  in  reference  to  Eq.  (15)]  whether  or  not  the  small¬ 
est  scale  of  the  physical  process  has  been  resolved.  It  should  also 
tend  to  the  corresponding  limit  at  large  X  if  the  record  length  L 
spans  a  sufficient  number  of  outer  flow  scales.  Accordingly,  we 
can  always  expect  the  limits  in  Eq.  (16)  to  be  attained  by  the 
analysis  of  proper  data. 

Of  course,  the  preceding  alone  does  not  preclude,  a  priori,  the 
existence  of  a  D,(X)  that  osculates  a  constant  value  over  a  range 
of  scales  X. 


<  *m,n  <  *  <  <  8. 

in  the  case  of  isoscalar  surfaces.  That  is  rendered  unlikely  by  the 
additional  requirement  for  a  characteristic  scale  X„  appropriate  to 
that  range,  as  noted  above. 

Rewriting  Eq.  (2).  we  conclude  that,  in  general,  we  must  allow 
the  exponent  D  to  be  a  function  of  X.  i.e.. 


Drf(X)  =  - 


d  log  Nd  (X) 
d  log  X 


(17) 


as  suggested  in  Miller  and  Dimotakis.  In  the  analysis  of  their 
turbulent  jet  fluid  scalar  data,  they  find  a  D,(X)  smoothly  rising 
from  0  to  1 ,  over  a  range  of  jet  Reynolds  numbers.  3.0  x  10  ^ 
Re  s  2.4  x  104.  with  no  indication  of  any  preferred  intermediate 
value.  They  also  report  that  the  observed  D,(X)  is  very  close  to 
one  derived  from  a  simple  model  process  in  which  scalar  thresh¬ 
old  transitions  are  log-normally  spaced. 

Figure  1  depicts  TV, (X),  as  required  to  cover  transitions  of  the 
scalar  jet  fluid  concentration  trace  c(t).  of  a  threshold  set  at  the 
local  mean  c*  on  the  axis  of  a  turbulent  jet  in  water,  at  x/d  = 
100,  and  a  Reynolds  number  of  1.15  x  !04.  In  this  plot,  the 
covering  element  scale  X  was  converted  to  an  equivalent  spatial 
extent,  measured  in  meters,  to  facilitate  comparison  with  other 
temporal  data,  at  both  higher  and  lower  Reynolds  numbers  (cf. 
Miller  and  Dimotakis  [1991,  Fig.  4]).  The  outer  scale  (local  jet 
diameter)  is  estimated  to  be  8  —  0. II  m,  i.e..  Iog,0(8/m)  — 


'Note  reference  to  the  use  of  Taylor's  frozen  floss  hypothesis  in  Sreenivasan. - 

Ramshankar  and  Meneveau.  ( 1 989.  Table  l|  Note  Id),  where  the  “The  important  issue  of  the  choice  of  the  threshold  was  addressed  in  Miller  and 
Sreenivasan  and  Meneveau  1 1986)  jet  data  are  cited  Dimotakis  1 1991)  The  reader  is  referred  to  that  discussion  for  further  details 
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-  1.0.  As  can  be  seen,  the  computed  coverage  observes  the 
required  limits  noted  in  Eq  (15). 

Figure  2  plots  the  logarithmic  derivative  D,(\).  computed  from 
the  /V,(\)  coverage  plotted  in  Figure  1.  Note  that  the  measured 
D,(X)  — »  0.  at  the  smallest  scales,  and  D,(A)  — »  1  for  A  ~  8.  or 
larger,  as  required.  The  dashed  line  in  Figure  2  marks  the  level 
D,  =  0.36.  The  latter  was  the  value  reported  to  have  been  de¬ 
rived  from  direct  measurements  of  one-dimensional  temporal 
data  by  Sreenivasan  and  Menevcau  1 19861.  It  would  also  be  the 
value  derived  from  the  previously  cited  estimate  of  D,  =>  2.36  if 
the  assertion  that  O,  =  D,  +  2  is  accepted  |Eq.  (I3)|.  However, 
while  the  assertion  of  Eq  ( 13)  may  be  true  in  special  cases,  it 
must  in  general  be  false,  as  will  be  argued  below  . 

Consider  an  isotropic  random  surface  in  three  dimensions, 
which  we  accept  as  possessing  a  fractal  dimension  D,.  and  its 
intersection  with  a  plane,  yielding  a  two-dimensional  curve, 
which  we  will  also  accept  as  possessing  a  fractal  dimension  D: 
Using  the  two-dimensional  curve,  wc  generate  (extrude)  a  cylin¬ 
drical  surface  by  translating  the  two-dimensional  curve  along  a 
direction  perpendicular  to  its  plane  As  Fcder  |I988,  Fig.  13.1 
and  related  discussion)  notes,  the  dimension  of  the  resulting  cy¬ 
lindrical  surface,  which  has  very  different  space-filling  proper¬ 
ties,  is  given  by  D2  +  I  Equivalently,  we  can  expect  the 
required  coverage  fraction  F,(A).  cf.  Eqs.  6.  8.  and  9.  to  be 
different  for  the  original  three-dimensional  isotropic  surface,  vs. 
the  generated  cylindrical  surface 

This  suggests  that  for  a  random  surface  in  three-dimensions,  of 
dimension  D,.  and  a  planar  section  of  the  same  surface  of  di¬ 
mension  D2.  we  can  expect,  in  general. 

Dy  *  D:  +  I  (18) 

(analogously,  for  curves  that  are  originally  planar  of  dimension 
D:  versus  linear  cuts  yielding  intersections  of  dimension  O,). 

Is  the  implication  that  there  can  be  no  (power  law)  fractals? 
Not  at  all. 

To  return  to  the  earlier  point  on  similarity  an  mensional 
analysis,  the  inference  that  should  be  drawn  is  that  a  (power  law) 
fractal  coverage  relation  requires  a  characteristic  scale.  Without 
that,  and  as  opposed  to  the  original  proposals,  we  cannot  write 
Eq.  (1)  dimensionally  correctly.  Conversely,  processes  that  do 
possess  relevant  characteristic  scales  could  exhibit  (power  law) 
fractal  behavior. 

What  is  then  the  relevant  characteristic  scale  for  the  triadic 
Cantor  set  (e.g..  Mandelbrot  |1977,  p.  98|).  often  described  as 


self-similar,  for  example? 

It  is  easy  to  appreciate  that,  for  the  triadic  Cantor  set.  the 
characteristic  scale  is  the  original  length  L  that  is  successively 
subdivided  into  thirds,  i.e..  \0  =  L  in  Eq.  9.  The  resulting  set  is 
obviously  not  self-similar  at  ail  scalings:  it  is  not  the  same  if  we 
cut  the  original  length  L  in  half,  for  example.  It  is  also  not 
self-similar  homogeneously:  successive  subdivisions  do  not  en¬ 
dow  the  middle  third  of  the  original  L  with  any  Cantor  "dust." 
Only  the  first  and  last  third  subsegments  of  any  (sub (segment  of 
the  Cantor  set  are  similar.  This  represents  a  rather  restricted 
similarity  scaling,  which  is  pegged  to  the  original  length  L  and  a 
coverage  phased  with  respect  to  the  original  segment.  The  oscil¬ 
lations  in  the  calculations  of  Smith  et  al.  (1986).  as  well  as  in  the 
spectrum  of  the  triadic  Cantor  set.?  can  be  understood  in  this 
light. 

Accepting  that  the  geometry  of  coastlines  may  be  describablc 
by  a  power  law  fractal  relation,  for  example  (cf.  Feder  1 1988. 
Fig.  2.7)).  one  can  think  of  several  candidate  characteristic  scales 
that  could  be  relevant.  Potentially,  these  range  from  buoyant 
scales  of  the  motion  w’ithin  the  earth's  mantle,  responsible  for 
tectonic  motion  on  continental  dimensions,  to  ocean  wave  scales 
at  intermediate  topography  scales,  to  rock  grain  boundary  and 
fracture  scales,  and  below,  as  appropriate  for  microscopic  and 
submicroscopic  stochastic  geometry  scales.  One  could  imagine 
the  application  of  Eq.  (9).  corresponding  to  each  of  these  separate 
characteristic  scales,  with,  potentially,  a  different  value  of  the 
exponent  D  for  each  one. 

In  the  context  of  fluid  mechanics  and  turbulence,  this  is  also 
not  to  say  that  (power  law  )  fractal  relations  and  constant  fractal 
dimension  exponents  are  ruled  out.  The  arguments  here  attempt 
to  address  issues  arising  in  fully  developed,  turbulent  flow  in  the 
intermediate,  inertial  range  of  scales,  where  power  law  fractals 
have  been  proposed  to  apply.  Conversely,  however,  isoscalar 
surfaces  that,  for  example,  demark  the  boundaries  of  a  turbulent 
region  and  even  may  serve  as  the  definition  of  the  outer  scale  b 
are  susceptible  to  dynamics  that  arc  obviously  not  independent  of 
a  characteristic  scale,  namely,  the  outer  scale  5.  Such  interfaces 
would  then  appear  to  possess  the  necessary  condition  for  a  power 
law  fractal  description  to  be  applicable  and  may  indeed  be  found 
to  display  such  behavior.  It  would  then,  however,  be  surprising 
if  the  resulting  fractal  dimension  exponents  were  found  to  be 
universal,  i.e..  independent  of  the  details  of  the  How  geometrv. 
Reynolds  number.  Schmidt  number,  etc. 

More  generally,  however,  one  might  accept  the  need  to  ana¬ 
lyze  stochastic  geometric  data,  using  the  tools  of  fractal  analysis, 
without  requiring  that  they  conform  to  power  law  similarity  re¬ 
lations.  Non-power-law  stochastic  geometry  similarity  could  be 
useful  for  analogous  purposes,  as  was  suggested  by  Miller  and 
Dimotakis.  who  instead  found  log-normal  similarity  applicable  in 
their  analysis  of  turbulent  jet  fluid  scalar  fluctuation  records. 

Indeed,  the  log-normal  distribution  was  originally  proposed  by 
Kolmogorov  (1941b)  to  describe  the  self-similar  fragmentation 
process  resulting  from  successive  rock  crushings. 
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